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- ZEITSCHRIFT FÜR ANGEWANDTE 
MATHEMATIK UND MECHANIK 
INGENIEURWISSENSCHAFTLICHE FORSCHUNGSARBEITEN 


'HAUPTAUFSÄTZE 


Verzerrungstensor, Verzerrungsdeviator und Spannungstensor 
bei endlichen Formänderungen 
Von Hans Richter in Haltingen/Lörrach 


Ps werden Postulate aufgestellt, denen bei der Bildung des Verzerrungstensors, des Verzerrungsdeviators 
und des Spannungstensors zu genügen ist, und hieraus die allgemeine Gestalt dieser Tensoren in beliebigen 
Koordinaten abgeleitet. Als einfachste Definition des Verzerrungstensors erscheint die gemischt-variante 
logarithmische Deformationsmatriz, wo der Deviator in üblicher Weise gebildet werden kann, und wo die In- 
varianten des letzteren die Beanspruchung invariant charakterisieren. Bei entsprechender Definition des Span- 
nungstensors bleibt die Gestalt des allgemeinen Elastizitätsgesetzes invariant gegen Koordinatentransformation. 


Postulates are laid down that have to be satisfied on forming the distortion tensor, the distortion. deviator, 
and the tension tensor, and thus the general form of these tensors are deduced in arbitrary co-ordinates. T’he 
combined-variant logarithmical deformation matrix proves the simplest definition of the distortion tensor. T’he 
deviator mar be formed in the usual manner, and the invariants of it characterize the strain in an invariant 
way. If the tension tensor is defined accordingly, the form of the general law of elasticity continues to be in- 
variant to co-ordinate transformations. ; : 


On etablit des postulats pour la formation du tenseur de deformation, du deviateur de deformation etdu 
tenseur de tension. La forme generale de ces tenseurs en coordonnees arbitraüres en est deduite. La matrice lo- 
garithmique (mixte-variante) de dejormation fournit la plus simple definition du tenseur de deformation. 
Le deviateur peut &tre fjorme comme de coutume et ses invariantes caracterisent la sollicitation. d’une maniere 
invariante. Le tenseur de tension &ant defini conformement, la forme de la loi generale d’elasticite reste 
invariante dans toute transformation .je coordonnees. 


VEeTAaHaBAUBAOTCA TIOCTYIATEI, U3 KOTOPEIX HY3KHO UCXOANUTB IHPH NOCTPoeHHH TeH30Pp&a 
zehopmaunf, neBuaropa nebopmannı u TeH30pa HANpMKeHuH. Orc onpenensetca o0maa 
dopMa 9TUX TEeH30POB B IIPOHSBONBHEIX KOOPAuHaTax. Hanöonee NpocTsim onpere,renneM 
TeH30pa NehopMmannä  MpencTaBgugeTcH CMEINAHHO-BAPHAHTHAA NOTAPHPMUYECKAA. MATPuLa 
zebopmanuf, MpmyeMm ANeBHaTop CTPOHTCA OÖBIUHBIM c10OCO60oM. MHBapmaHTsı NeBuaTopa 
HHBAPHAHTHO XapakTepusymT HaupsweHusa. Ilpıu COOTBeTCTBEHHOM OLpeneleHuN TeH30PA 
HAIPASKeHmÜu dOPMA OCHOBHLIX YPpaBHeHnÄ TeOPAHM YIPYTocTu MHBAPHAHTHA IIO OTHOMIEHHIO 
.E Hpeo6pa30BaHnIo KOOPAHHAT. 

$ 1. Einleitung 


In der Theorie endlicher elastischer oder plastischer Verformungen geht man im allge- 
meinen von demjenigen Verzerrungstensor aus, der für allgemeine Koordinaten durch Bildung 
der Differenz der Quadrate der Linienelemente im deformierten und im Ausgangszustand ent- 
steht 2). Die Verwendung gerade dieser Charakterisierung des Verzerrungszustandes ist natürlich 
nicht zwingend vorgeschrieben. Im Gegenteil erscheint bei einer genaueren Analyse, die im 
folgenden durchgeführt werden soll, gerade diese übliche Definition des Verzerrungstensors 
nicht als diejenige, die dem Problem der Untersuchung endlicher Verzerrungen besonders gut 
angepaßt ist. So führt bereits die Aufgabe, aus dem üblichen Verzerrungstensor einen Deviator: 
abzuleiten, der in Absonderung der Volumänderung nur die Gestaltänderung charakterisiert, 
zu .eigenartigen Schwierigkeiten und Mehrdeutigkeiten !). Der tiefere Grund hierfür liegt wohl 
darin, daß man sich bei der Behandlung der endlichen Formänderungen zu stark an das Vorbild 
der infinitesimalen Verzerrungen gehalten hat, bei denen man eine beliebige Deformation durch 
additive Aufspaltung in den symmetrischen und den schiefsymmetrischen Anteil in eine reine 
Streckung und eine reine Drehung zerlegen kann. Bei endlichen Verformungen ist diese additive 
Aufspaltung jedoch nicht mehr möglich; an ihre Stelle tritt eine multiplikative Zerlegung der 
allgemeinen Deformation in eine Drehung und eine Streckung, wobei diese Faktoren nicht mehr 
kommutativ sind. Jeder Versuch, Definitionen durch additive Zerspaltung zu bilden, muß 
daher auf grundsätzliche Schwierigkeiten stoßen. 

Wir wollen nun in dieser Arbeit so vorgehen, daß wir — gewissermaßen axiomatisch — 
an die zu bildenden Definitionen von vornherein gewisse Forderungen stellen, die uns als zweck- 
mäßig erscheinen, und dann zeigen, wie aus diesen Forderungen sich bestimmte Möglichkeiten 
für diese Definitionen als besonders naheliegend ergeben. 


1) Vol. R.Moufang: Volumtreue Verzerrungen bei endlichen Formänderungen. Z. angew. Math. 
Mech., Bd. 25/27. (1947), 5. 209—214. 
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zu 4. 


25 Es N 6 af} = [Bf (MW) aa), falls BA—= AB, jedoch nicht notwendig B-dA 


.“ ER " x 

1. Mit großen gotischen Buchstaben 4, ®, ... bezeichnen wir « hige 
an 2), 0 (Mir ist das Element in der ö-ten Zeile und der k-ten. Era a 
minante von X. {X} ist die Spur von %, also die Summe der Elemente der F PaR 
ist die an der Hauptdiagonale Gespiegelte von WU. € ist die Einheitsmatrix. ist die In 


2. Kleine gotische Buchstaben t, 4, ... bedeuten Vektoren: i= (ar er 2); u 
xt) ist das innere Produkt. x xy ist das äußere Produkt. 


3. Ar entsteht durch Anwendung von W auf x: (Ar) = Zaun 2 


4, Produkte BA sind von rechts nach links zu lesen: (BWL— B(Ar). 
5. Ist f(a) = N db, x", so ist unter Voraussetzung der Konvergenz: f (A) — — EM; ER 
ara) = FU HIN) — 0, was mit F(WAA nur bei W-dY=dX- A übereinstimmt 9.77 


b) Hilfssätze 
(2.1) {U,Uz... An} bleibt bei zyklischer Vertauschung der Faktoren ungeändert. 


(2.2) Alle Invarianten von X gegen afline Transformation > EAC- sind Funktionen’ der drei 
Invarianten j= {[W, k= {W} und, 7= {W}. Die charakteristische Gleichung von 4 ist: 


; Urs 1 1 
Bath SjR44P)=0. 


(2.3) "Es ist . FCACH = EM CE. 
(2.4) Hat X positiv reelle Eigenwerte, so ist In X definiert. Es ist {In W} = In |. 


—=dA-Bi 


(2.6) In a Koordinaten ist eine reine Sue © symmetrisch mit positiven Eigen- 
werten. 


(2.7) In kartesiechen Koordinaten gilt für eine euklidische Transformation R: AR = €. 


(2.8) Ein beliebiges A mit |X| #0 läßt sich eindeutig darstellen in der Form Y= &- N, d.h. 
als euklidische Transformation mit nachfolgender reiner Streekung. Bei U] > 0 ist R eine 
direkte Transformation, also eine reine euklidische Drehung. 


(2.9) Es ist 2A =dH-Ar. j 
(2.10) Es sei y = Ux eine Koordinatentransformation, bei der X in A* übergeht. W ist ein 
zweifach kontravarianter Tensor, falls A* = UAU- Ju], 


zweifach kovarianter Tensor, falls A* = NIAU- u", 

kontravariant-kovarianter Tensor, falls A* = UAU- |uj®, ; 

kovariant-kontravarianter Tensor, falls A* = UFAU- [Ul® ist. 
j 


Bein=0ist Wein eigentlicher Tensor; bein # 0 eine tensorielle Dichte. (Die Überein- 
stimmung dieser etwas weniger geläufigen Darstellung der Tensoreigenschaft mit der üblichen 


ergibt sich unmittelbar, wenn man symbolisch (N); = x: y; setzt, wo x und t) kontravariante » 
oder kovariante Vektoren sind). 


(2.11) Essei X = Urund y’ = Uy; dann ist! xy = [U UF xy). 


$ 3. Der Verzerrungstensor 


Wir wollen uns nun überlegen, welche Forderungen wir billigerweise an die Verzerrungs- 
matrix stellen können, um dann die Realisierbarkeit dieser Forderungen zu untersuchen. 


Es sei X die Matrix, die die Umgebung eines Punktes £ auf die Umgebung des Bildpunktes x 
abbildet: 


(3.1) dr= Ad}. 


?) Ob eine Matrix ein Tensor ist, ist an (2.10) zu sehen, 
») Vgl. jedoch (2.5). 


| t ‚den erreichten a an. Bei Se Material, wo der Spankta 
zustand nieht nur vom erreichten Verzerrungszustand abhängt, sondern auch von dem Weg« 
Fe der zu diesem führte, ist, dann die Angabe von X allein nicht ausreichend. Bei elastiscom 
® N rial dagegen genügt A zur Charakterisierung. der Verzerrung. Bei anisotropem Material. FE 
7 ist auch die in W enthaltene. Drehung wesentlich. Es ist dann A selbst zur Beschreibung der 
' ER Verzerrung zu verwenden, während alle Verzerrungstensoren, die wie der gewöhnliche eine 

: "euklidische Drehung eliminieren, unbrauchbar sind. Die Aufstellung: von alten Ne IS en 
" Ban hat daher überhaupt a für ‚isotropes. Material Sim. cv Pa 


| REN x ne RENNEN THE RR a) Dostulate R I ee BR ER 
EN D mgemäß soll nun unter der ausdrücklichen re der Anwäntbarneil 
isotropes Material zu X eine Verzerrungsmatrix ® (X) definiert werden‘). Während nun 
sicher kein Tensor ist, da X auf zwei verschiedene Punkte bezogen ist, wollen ‚wir die Tensor- 24 
. eigenschaft von % fordern. ‚Damit ergibt sich das erste Postulat: | Br \ 


SER E Bist ein Tensor, der aus W und den Matrizen der Metrik TI 
a E undrgebildet werdenkann. | 


i Weit rhin soll bei 8 die. unwesentliche i in en Drehung unberneksichtigh bleiben. 
"Es ‚soll sich also % nicht ändern, wenn vo r "Ausübung von X erst eine euklidische Drehung R_ = 
in 2 durchgeführt wird. Man kann sich statt dessen auch auf den Standpunkt stellen, daßeine 
nach Xin x durchgeführte Drehung den Verzerrungstensor nicht beeinflussen soll. Dies würde 
AR bedeuten, daß X als Verzerrung | in ? mit nachfolgender unwesentlicher Drehung aufgefaßt wird. 


| Den so zu t gehörigen Verzerrungstensor wollen wir mit ® bezeichnen. Die Analyse von ® ‘ : 
‚und 3 R ist vollkommen analog, so daß wir uns im folgenden auf die von % beschränken und die ie 


E : für ® analogen Ergebnisse lediglich vermerken, wobei die entsprechenden Größen durch ein * 
= gekennzeichnet werden. 


Be EENN, Die genannte Eigenschäft v von ® und ® drückt sieh nun aus durch das Bosalalt 

v2: VAR) = BM, resp. DR A) = RW. 

ee Weiterhin verlangen wir noch ein Superpositionsprinzip für koaxiale reine 
7 Streckungen durch Postulat 


E25 #3: EstseieH S, und ©, zwei Buakiela Streekungen: = 58. 
Er Es sei %, = % (6), = %(5) und B=R(6,6,. Dann soll es eine 
A ee umkehrbar eindeutige Funktion /f(») geben,.so daß fı(B) 
+f(B,) = f(B) ist. f kann gegebenenfalls abhängig vom Koordinatensystem sein. al 
Schließlich müssen wir noch verlangen, daß für infinitesimale Verzerrungen die neue - 
Definition in die alte übergeht. Dies liefert die Limesbeziehung 


v4: Für infinitesimale Deformatio nen“ E+4B a kartesischen ü ee 


B- Koordinaten Scht.der TEN SER UNE in zAB+AS) Fo (8) 
über °). | 


{ 
/ 
\ 


b) Die Realisierung der Postulate in kartesischen Koordinaten 


En Aus Vereinfachungsgründen wollen wir zunächst kartesische Koordinaten voraussetzen. 
2 _ Original- und Bildpunkt der Deformation nennen wir dann $ und y. Die Deformationsmatrix k 


- heiße jetzt ®. Die zugehörigen Verzerrungstensoren sind W und ®. 
Gemäß (2.8) schreiben wir zunächst 
(3.3) | B-SR-NE bi G=-HISN. 
Zur Auffindung dieser Zerlegung bildet man zunächst ® 8. Nunist fürt£0:0< Br: 8% 
was mit Hilfe von (2.9) liefert: 0<r'® ®r. Die symmetrische Matrix ® ® ist also positiv IR 
| E 4) Die Bezeichnung ® (X) soll dabei nicht heißen, daß ® eine Funktion von WX im Sinne von $ 2. Pkt.5 i 
ist, sondern nur, daß ® zu X gehört. x: 


q 
5) Wie üblich bedeutet y=0(2): lim = = 0. 


al A Dr ee 


 efinkk And hat Aalır erndenig: fe peitie definite ( 
‚sich dann als R= &1%. Entsprechend ist &—! 2. AR 


.  tensoren beschränken, die zu.reinen Streckungen gehören. BR er IR | 


=d&+0(d&) und ®(E+1dC)=4d&+oLdS), we 


3=f(W), so wird damit für infinitesimale Streckungen: 3 (€ +46) 


' Eigenwerte von © bilden: , es 


Es ist dann: 8-1} und damit für große n: Ve=C+l2+0(t ji Er, Be | 
sye)- ne ). Nach'-7.:8 st’ welter BIETE ar A U 8 +n:0(4). Da = 


| (3,5%) 012), x +0(e); Ari+)=xz+o0(ed); kil 14924, +o(e). 


Be.‘ 
L o 


Nach V2ist W(B) = W(S), resp. Bd) = BIO). Wir können u sal ie, 
'Es seinun & eine infinitesimale Streckung: ehe Day ist mach Pa: | Re R 
Ede del Sr 
en &) u | x 
Sei nun wieder & eine endliche reine Streckung, so läßt sich N ge post 


Aus 7 3 ergibt sich dann f(d& +0(4©&)) +f(AdS +o(d Die 
Da dies für jedes A und d& gelten muß, folgt für kleine «: f(«)= x +0 (z) 


(3.4) 8—1n ©; ae 8 In &: „logarithmische Delormakiannmakrheng 


die linke Seite dieser Gleichung von n unabhängig ist, können wir n gegen unendlich gehen lassen ee 
und erhalten: 8(©&) = 2. Hieraus folgt insbesondere, daß f(x) bis auf einen willkürlichen 
Faktor eindeutig bestimmt ist. 


Ist nun 9 die Umkehrfunkton zu f, so haben wir schließlich, da ja 2 eine umkehrbar ein- 
deutige Funktion von © und damit auch von &=B%8 ist: 
35) BEI H-AHOS-EBB; resp. B-IG-Kl(- BR). 


Der Ansatz (3.5) genügt umgekehrt stets den Postulaten Y2 und V 3, wobei f eindeutig als 
Umkehrfunktion von g zu wählen ist, während die Erfüllung der Limesbedingung V 4 fordert, 
daß für kleine x gilt: 


In der Tat ist dann für infinitesimale Verzerrungen B= & +48: 
BÜ-ELHABHLE) HB),  S-VBÜ-E+,ABLAB)+0(AR), 
&=In ©- 548 +AB) +0(dB). 


Es ist also bei jedem (3.58) genügenden g: ® = = (dB +4) +o(dB). 


Jeder mit unseren Bee verträgliche Verzerrungs- 
tensorist damitals Funktion derlogarithmischen Deformations- 
matrix erkannt. Vom Standpunkte unserer Postulate aus erscheint ®W—= & als die 
einfachste Definition der Verzerrungsmatrix, da hier das Superpositionsprinzip mit f(x) = x 
erfüllt ist. Wir werden später sehen, daß auch für die Bildung des Deviators diese Defini- 


tion als die einfachste erscheinen wird. 
Wird anschließend an X noch eine euklidische Drehung N, durchgeführt, so geht wegen 
RAFENKSR-N, ORT N,R die Streckung & in R, SH- über. Den gleichen Über- 2 


gang haben wir, wenn wir eine euklidische Koordinatentransformation y,= R,y durchführen. 
Gemäß (2.3) geht dann Win A (RN, SORT) = N,WNZ! über. Die Achsen von ® werden also 


. bei nachträglicher Anwendung von N, einfach mitgedreht. Bei der anderen Auffassung der 


letzten Formel als Ergebnis einer Koordinatentransformation sehen wir an (2. 10), daß ®W sich 


wie ein Tensor transformiert, wobei wegen NR! — R, kein Unterschied bezüglich der Varianz 
vorhanden ist. Entsprechendes gilt natürlich für R. 


c) Erweiterungaufkrummlinige Koordinaten 
Wir gehen jetzt von den kartesischen Koordinaten 4 auf beliebige Koordinaten x über: 
2?=2(y). Für eine Umgebung des undeformierten Materials seid? — üad, für die entsprechende 
Umgebung im deformierten Material sei dt = Udy. Ü und U sind die Funktionalmatrizen 
vonTt=r(h) DD p.nnon ’ 


°) Da mit f jedes Vielfache von f dem Postulat 7 3 genügt, kann man fo) zu 1 normieren. 


_ 


gegangen in. 
Bra... RE UR, Mu “ 
2 Umgekehrt ist RN 
| De a sam 
B de ‚wir sofort erhalten: R 
RE a gehe ums, 
Br | Mar &=BBd= nuasal. 


E 
= 
4 
2 
E 
. 


N RG, Sin Tohtenelene 
ae ae Ai REN ER ER 


le Matrix der Metrik in Entsprechend de dee N BEREL 
EB Se a ET 
die Metrik a ' en | ER 


wir mit, Ztte von 2 29: m 


hr „ 
NH, 


Die Deformation des Materials erscheint jetzt als: 4 = -u8 im a. ® ist also übers. 


Mit Hilfe der beiden. fetten Formeln können wir also die zu den SE ee Koordinaten A 


. gehörigen Matrizen 8, S und & durch X und die Übergangsmatrizen U und ü ausdrücken. 


a) Falldes ingewischten Tensors 
Wir nehmen zunächst an, es sei der Verzerrungstensor ® zweifach kontravariant definiert 
und genüge den Postulaten V 1 bis V4. Dann muß nach (2.10) sein: 


V-UBU, 


| wobei % einer der Teroren (3.5) ist. 


Um die besondere Gestalt von W kennenzulernen, nehmen wir des Spezialfall, daß 2 
eine reine Streckung © in den Koordinatenachsen und U ame koaxial ist, also \ 


.&=| 4 guda:.Me= 005 } Wr 


Dann ist wegen (3.5) ® ebenfalls auf Hauptachsen mit den Eigenwerten 0, h(A,). Das Super- 
positionsprinzip fordert nun die Existenz einer Funktion f(x), deren Koeffizienten die g, ent- 5 
halten können, so daß gilt: 


(3.9) | SR (a)) + hm) =F lei h (Am) 


für beliebige A, und u,. Es ist. daher 


(ha) = C,' "Ind, 0,= 0,(01 0» 03)» 


Durch Differentiation nach A erhalten wir hieraus: 
(3.10) (ha) WM) 0,7: 


Setzen wir speziell A=1, ‚so wird wegen (3.5%): o,'f’ (0) = (,, so daß bei der Normierung 


FE Ws = 1 wird: ; 


f(&h (A) ni 1-W ()' 


ae 
Die rechte Seite dieser Gleichung ist unabhängig von den e,. Es muß daher f(x) =f’(0)=1 


und damit A (A)=InA sein. Es ist dann ® = In S= % und damit 
Bg=ULU 


zu setzen. Umgekehrt erfüllt diese Definition für ® offenbar auch bei beliebiger Wahl von 8 
und ll alle Postulate V,1 bis V 4, da ja für 2 das Superpositionsprinzip rein additiv ist und sich 


‚daher bei Multiplikation von links mit U und von rechts mit U einfach auf ein additives Gesetz 


für ® überträgt. . 


eine mit unsgren Postulaten BER deti 
keit: B- yalık 


SR N U-- um ee iu. 
| Aline ergibt sich“ 


N ‚in der Tat ® und ® nur von X, Ö und & abhängen. Die Darstellung durch diese Matrizen ist 
jedoch sehr umständlich. Außerdem sehen wir, daß die Invarianten von ® andere sind als. re 
" von ®. Dies ist insofern unangenehm, als z. B. in der Elastizitätstheorie endlicher V: 


- mation ändern, so ergeben sich hieraus zusätzliche ‘Schwierigkeiten. Auch ist es dann nicht Krach 2 


| Wegen (2.3) genügt dann ® automatisch dem Superpositionsprinzip mit dem gleichen f(x) 


Es ER US6 für Te, en soriel 


$ m ano. | E 
Entwickeln wir den In für nicht zu große Verzerrungen in eine Potönzibihe, so sehen an a I 


anzunehmen ist, daß die thermodynamischen Größen wie innere Energie, Entropie usw. Funk- 
tionen der Invarianten der Verzerrung sind. Wenn dieselben sich nun bei Koordinatentransfor- 


mehr möglich, die Beanspruchung mit Hilfe der Invarianten BunphBHDE von der Kr 


wahl zu beschreiben (vgl. $ 4). ’ Fan | 


Die entsprechenden Betrachtungen gelten natürlich für den Fall, daß ® oder ® zweifach‘ a 
kovariant ist. Es wird daher zweckmäßig sein, auf den mit den ungemischten Tensoren ver 
bundenen Vorteil, der Symmetrie zu verzichten. 


ß) Falldes gemischten Tensörs Si 5 
In dem Falle, daß ® kovariant-kontravariant ist, muß nach (2.10) sein: B- UIWU. 


« 
kenne 


wie ®. Insbesondere ist also das eindeutig bestimmte, normierte f(x) unabhängig von der 
Koordinatenwahl. Weiterhin hat ® dieselben Invarianten wie W. Jede Funklion von ®, deren 
Koeffizienten von den Invarianten von ® abhängig sind, überträgt sich auf die gleiche Funk- 
tion von ®. 

Aus der einfachsten Definition W— 2 erhalten wir jetzt bei beliebigen Koordinaten: 
2* = 1121 oder wegen (3.4) und (3,8): 


- . 
Tg a u wu nr 


1 — 2 an 

t=7 An UIWNETAUHU 
oder # 
1 Al | 
(3.12) tin (SEAGTIA): „logarithmischer Verzerrungstensor“, i 
2* gehorcht dem Superpositionsprinzip mit f(a)=x. ö i 
1 
Der allgemeinste unseren Postulatengenügende Verzerrunge- AN 
tensoristdann: ee | 
(3.13) | 89 (2%) 7 EC DER ENG Ä 
wog, hundkden Bedingungen (. 5%) unterliege k 


Ganz analog erhält man 


— 
A 


” 


8-9) = YAHAS) -RÜHUSH, = Im ASAG-). 


Bis auf einen beliebigen Faktor ist die Funktion /(x) des Superpositionsprinzips die Umkehr- 
funktion von = g (Y). 


Ganz entsprechend ist vorzugehen, wenn ® und ® ‚kontravariant-kovariant sein sollen. 
Es wird dann 


Ur RR 


gr = z In (AA 6) 
(3.12%) und | 
| gr — x n(GIAHW. 


Alle anderen Beziehungen bleiben ungeändert. 


TEEN 


‚und entsprechend 


N wi Kern dh na EN 0 a en I 


- 
t 


Richter, Verzerrungstensor, Verzerrungsdeviator und Spannungstensor 71 
< 5 ; \ ü 


d) Berechnung der Volumdehnung» 
Die mit WA verbundene Volumdehnung ist v = ||; also nach (3.7*): 
| v— uf [ll al. 
Andererseils ist nach (3.6), (3,12) und (3.12%): 
jer|j=|®]: SH - [Up = 22. 
Also ist nach (2.4): 
(3.142) na = der ; 
oder nach (3.15) 


(3.14) , e= {fi} = FR 


e) Beziehung zumgewöhnlichen Verzerrungstensor 
Die gewöhnliche Definition des Verzerrungstensors T, resp. €, ist: bekanntlich Ai 
de —d2—=2:dr-Tdr=2dt- Taf. 
Nun erhalten wir mit Hilfe von (2.9): 
d2=dı: Gdr=AdE- GAdFT=AdLE- AGAIE, 


d2—= dy- ATGAIdr —=d}-Gd?}. 
Also wird: 

2T=- B— (AGAIN und 2T=-AGA—G. 
Um die Varianz zu sehen, formen wir nach (3.8) um: 


(3.15) 2T=-UL(E- UATEAIU)- Ut= UF (E69 -U- 
und entsprechend x 
(3.153) 2-12 &@ _- Hü-. 


Nach (2.10) sind T und T also zweifach kovariante symmetrische Tensoren. Das Superposilions- 
prinzip ist für sie nicht erfüllt und die Invarianten ändern sich bei Koordinatentransformation. 


Dagegen genügen die gemischten Tensoren TO, 61T, 61T und T 6 allen gestellten 


1 
Postulaten Y 1 bis V 4. Für das Superpositionsprinzip ist nach (3.15) f(x) = — 5 In (1—2x) 
für T 6-1 und. GT, dagegen f(x) -; In(1 +2) für GT und TG- zu setzen. 


Für die Volumdehnung » ist dann nach (3.14): 
2 = |E-2T641- [E-26274= E+2 628] = |E +2T67. 


$4. Der Verzerrungsdeviator 
a) Postulate 

Der Verzerrungsdeviator ® soll aus dem Verzerrungstensor abgeleitet werden und ledig- 
lich die mit der Verzerrung verbundene Gestaltänderung charakterisieren. Damit ergeben sich 
bereits die an ihn zu stellenden Postulate: 3 
Di: Unterscheidensich zwei Deformationennurdurcheine Ahn- 

lichkeitsstreckung, so haben sie das gleiche ®. 
D2: Ist mit der Deformation keine Volumänderung verbunden, 

soistd=P8. 


b) Die Realisierung der Postulate 
Eine Ähnlichkeitsstreckung im undeformierten oder im deformierten Zustand hat die 
Gestalt A@, A > 0, mit der Volumdehnung 4°. Setzen wir 


DIR es »U3 & my 1/3 D — vU3 & * Y, 9 
so ist nach D1: D(A) = D(A,). Da A, mit keiner Volumdehnung N ist Se 
DM) — B(U)—=9(8%), wobei nach (3.12) und (3.12%) gilt: = — „mv: & +%*. Hier- 


2) Vgl. 2. B. Moufang 1. ce. 


N ds önoen wir EDER 


gr 


‚Ganz ae ist 


| | a N = er 1 En 
esechrt sind bei dieser Definition offenbar auch die Re D 1 a D: ra. Mi 
nämlich U mit A > 0 multipliziert, so geht &* in &* +1nA- E über; &* und damit bleib | 
also ungeändert. Ist weiter »= 0, so ist nach (3.14°) {2 N also wog und d Pre , 
8 =g(2*%)=®%. Man beachte übrigens, daß ® automatisc Tensor ist, wenn wir Bals 
gemischten Tensor verwenden. Auch aus dieser Tatsache ereint Siehe eine BSR der 
gemischten Tensoren vor den ungemischten. n 
Am einfachsten wird die Deviatorbildung bei 8 = 2*, wo einfach ®= 2* wird. Die N 
Verwendung der logarithmischen Deformationsmatrix gestattet also auch bei BEUEDIGEN Ko. 
ordinaten die gewöhnliche Deviatorbildung. 
Es sei noch bemerkt, daß bei infinitesimalen Verzerrungen in kartesischen Kosriindien? 
‚der neue Deviatorbegriff in den alten übergeht. Ist nämlich B= € +d% eine infinitesimale 


Verzerrung, so ist 8=, (48 +48) +0(dB), so daß (4.2) wegen (3.5°) liefert: 
a ee 


x: - 1% mt 
a8 SON TE a ER +4B) +o(dB). 


re: 


I | 


Für den gewöhnlichen er Verzerrungstensor % eo” fanden wir in 33, e: Fin) 


Sa AR sh (1—22) also g(y)= — u —e2)) und v= JE -2T76722. Es it = 
E — m (E- 270) 6 - ne-rre- dur GE. 
FR Also wird schließlich ®): 
BE (4.3) d= |E-2TE1-W. [ze Cm an -1 - JE-2TB87]- e) 
Entsprechend ist 
(4.3) d= |E+2T 81-08 (f 4 — = Ne+2264)- 1- 6). 


c) Die Verzerrungsinvarianten 
Will man den Verzerrungszustand durch Invarianten charakterisieren, so wird man als 
erste geeignete Invariante von ® die Volumdehnung oder eine Funktion derselben wählen, _ 
während die beiden anderen Invarianten die Gestalländerung charakterisieren, also bei zu- 
sätzlicher Ähnlichkeitsstreckung ungeändert bleiben sollen. Da bei Benutzung der gemischten 
Tensoren — dies sei im folgenden stets vorausgesetzt — alle Invarianten von ® auch solche 
von %* sind, können wir nach (3.14%) als erste Invariante {%*} nehmen. Die beiden anderen 


Invarianten müssen nach Abschnitt (b) dann Invarianten von &* sein. Damit ergibt sich die 
Charakterisierung des Verzerrungszustandes durch 


| ee Eu! für die Volumänderung 
2 3 
y={%% ),2={8*} für die Gestaltänderung. 


Wegen &*— U EI ist auch y= {2%} und z= {&®}, so daß y und z die Gestaltänderung 
unabhängig von der Koordinatenwahl angeben, 


(4.4) 


e) Man vergleiche die etwas abweichende Bildung bei Moufang, 1, ce. 
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Wegen {f} = 0 heißt die Ey Base von & gemäß (2.2): 


4.5 _—<ı- 
| | ) x 5 % z 0. 
Damit diese Gleichung drei reelle Wurzeln hat, muß für 
ud 
(4.6) = ze } 
gelten: 
; 1 
(4.7) DR PER, IxI<z- 


Die geometrische Bedeutung von £ ergibt sich aus der folgenden Betrachtung. Es sei © eine 
beliebige reine Streckung. Dann können wir © auffassen als die n-malige Anwendung der 


« el 1 1 en 1 E7 
reinen Streckung ©, = ve. Dabei ist &, = In Y& Tr n6= Fr 2; also 2,= ru % und 
damit „= —,'yund„—=-,2. Hieraus ergibt sich (u = u =t. — Ist umgekehrt für 
zwei ce ©, und S: &, = 6, So ist y, = Ay und 2, =A°z. Gemäß (4.5) sind dann 
die Eigenwerte von ©, die A-te/Potenz der Eigenwerte von &,. Von einer Drehung abgesehen, 
ist damit ©, = ©. Wir können uns also ©, und ©, bis auf eine Änderung der Hauptachsen . 
als aus der gleichen infinitesimalen Streckung entstanden denken; bei negativem A mit Be- 
nutzung der Inversen. Dies bedeutet, daß £ die Beanspruchungs art charakterisiert. 

Speziell wird die einachsige volumtreue Verzerrung in geeignet gedrehten kartesischen 
Koordinaten dargestellt durch! 


SS UNE 3 
\ Sei FARO 
. 4 DI 
: 150 0 
Es ist dann t—2%t=1nA:[0 —1j2:°0 ; 
00. —12 
; 34 3 = 
also y=1n?}- > und: 21° r® womit sich ergibt: = re 
Dagegen wird für eine reine volumtreue Gleitung 
a nr 1+42 4 0 
1 er re) 
Bee Desert ) al 


Für die Eigenwerte von &?erhalten wir aus der charakteristischen Gleichung: A, ae Na a 
Für die Eigenwerte von gilt also: u, = 0, u; +us = 0. Damit wird y> 0,2=0; alsol=0. 
Wir haben somit gefunden: 
2 
[= charakterisiert die Beauspruchungsart, Der. eine 


ExtremwertZö=0 gehört zur reinen Gleitung und derandere Ex- 


tremwertc=- zur einachsigen Streckung. 


6 

Die Größe der Beanspruchung wird bei infinitesimalen Verzerrungen üblicherweise 
durch v{@% charakterisiert. Nun sahen wir, daß für infinitesimale Verzerrungen D = 8 ist; 
also gilt Yy als Maß für die Größe der Beanspruchung bei infinitesimalen Verzerrungen. 

- Ist nun wieder © eine endliche Streckung, so ist wie oben gezeigt: \y= N YYn- Da für 
genügend großes n yn die Größe der Beanspruchung bei Anwendung der infinitesimalen Ver- 
zerrung IS angibt, ist es vernünftig, Yy Zn Yun als Maß der Beanspruchungsgröße bei »-maliger 
Anwendung von VE. also für ©, zu verwenden. Damit haben wir: 


Yy gibt die Größe der Beanspruchungan. 


S5. Der Spannungstensor | 
Der Spannungstensor ® soll den Spannungszustand im Punkte x des deformierten Ma- 
terials beschreiben derart, daß bei passender Definition eines Flächenelementes d fin x die auf 


Spannungen direkt auf df zu beziehen, also ein ® zu konstruieren, ist physikalise 


 df wirkende Kraft durch‘ ®df gegeben wird. Wenn ı 
der Transformationsformeln natürlich auch in den Koordinaten von x ausdrüc 
“gang zu Lagrangeschen Koordinaten), so bleibt Pdoch an d gebunden. D 


Wir wollen daher darauf verzichten. _ ER 5 I, Se w Ak 
. a) Postulate Ce 


Bei kartesischen Koordinaten liefert die Spannungsmatrix ®, die auf ein Flächenelement Be 
df,im Punkte y wirkende Kraft 4f, in der Gestalt: df, = ®,d@f,. Im allgemeinen wird man | 


annehmen können, daß auf das Material keine solchen äußeren Kräfte wirken, die volumabhängige — $ 
Drehmomente hervorrufen. Es ist bekannt, daß dann ®, symmetrisch ist. Im folgenden braucht $ 


jedoch diese Symmetrie nicht vorausgesetzt zu werden. 


Haben wir beliebige krummlinige Koordinaten, so muß zunächst ein Flächenelement ar 


passend als Transformierte von df, definiert werden. Der Spannungstensor ® soll dann so 
gebildet werden, daß ®df wieder die Kraft angibt, die auf das Flächenelement ausgeübt wird. 
Wird das Element um den Vektor dz parallel verschoben ‚so soll die Arbeit da Pdf geleistet 
werden. Damit kommen wir zu den folgenden Postulaten: Y 
Pl: Bistein Tensorodereinetensorielle Dichte. 
P2: Für das Flächenelement gilt d{=@df,„ wo E noch passend zu 
bestimmenist. 


P3: Wird das Flächenelement um dz; verschoben, so wird die 
Arbeitd4A=dz- Pdfgeleistet. 


b) Die Realisierung.der Postlulate 

d A muß als numerische Größe invariant gegen Koordinatentransformation sein. Also gilt: 
(5-1) dz Bdiedze Bodfo> 
wenn dj, —= Udz der entsprechende Verschiebungsvektor in kartesischen Koordinaten ist. 
Nun erhalten wir mit Hilfe von P2 und (2.9): 

dzo: Bodfo= UNdz- PCrdi= dr EB UNdz—d;- UTB,C7gf. 
Der Vergleich mit (5.1) liefert, da d; und df beliebige Vektoren sind: j 
(3.2) R- UTC. 
(2.10) zeigt, daß P entweder &) zweifach kovariant beiC = U J ||" oder £) kovariant-kontra- 
variant bei € = U - Y|®|” ist. 

Im Falle (&) ist df kontravariant und im Falle () kovariant. Soll die Länge von df die 
geometrische Größe des Flächenelementes sein, so ist x= 0 zu selzen. ® ist dann ein echter 
‚Tensor; und zwar ist in den Fällen (a) und (P): 

(9,2%) R= u. RU 

und 
7} B-UTRT. 
Im Falle (&) ist dann df = Udf,. Wird das Flächenelement df, aus den Vektoren dy,o 
und dr,. gebildet, so ist df,=dr. Xdg„=Wtdr x Utdg, oder mit (2.11): 
df,=Y|®|-Uldr, x dxs). Damit wird 


(5.3°) di=Y|6|- 6 (dr, x dr). 
Im Falle (f) dagegen ist 
(5.3P) di= Y\®| (dr, x dx). 


Ob man mit kontravariantem oder kovariantem df rechnen will, ist natürlich gleichgültig. Wie 
(9.3) zeigt, liefert jedoch die kovariante Definition (8) die einfachere Formel. Allerdings ist 


dann nicht mehr symmetrisch, wenn ®, es war. Dafür bleiben jedoch alle Invarianten von B 
bei Koordinatentransformation ungeändert. 


c) Die Arbeitsleistung beiinfinitesimaler verzerruus / 

Wir nehmen jetzt an, es herrsche in einem Raumgebiete um ti) der homogene durch Bo 
definierte Spannungszustand. Ein geschlossenes Volumen V habe die Randfläche F mit den 
Flächenelementen df,. Wir führen nun in der Umgebung von t) eine homogene infinitesimale 
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Wieder "ergibt ‚die Ente yariänte Definition (B) die a Formel. 


-Der infinitesimalen: Verzerrung E+d%, entspricht in beliebigen Koordinaten die Ver- 
zerrung E+d% gemäß: U(E+4B,dy, = (E r d 2 Udy, also AB =UMTIBU 


. Wir haben daher mit (5.2) aus (5.4): 
; er ra ru nTanN - ussan). 
Be = m den Fallen [0 A (P) wird somit —_ ER 
4 9.4) RER AL Ga Bar) - u R 0 | ” £ 
4 Er = ER und f wre ( Er 
h ef) ae aA 1A Raw) = {a8 DE 
$ Wieder erhalten wir bei der Definition (ß) das einfachere Resultat. 
BE: d) Invarianz des Elastizität gesetzes 
nr "Für I Material hati in kartesischen Koordinaten das Elastizitätsgesetz die Gestalt 2 Re 

SE di= E42, 8 I Rrbeisieih), ke {0} Und {oo}, Be. 
RE, das elastische. Potential pro Volumeinheit des Ausgangszustandes 151.15 .4; | a. “= 
i nr Wollen wir erreichen, daß diese einfache Gestalt auch für beliebige Koordinaten gilt, so 
A _ müssen ® und 2 die gleiche gemischte Invarianz haben; denn nur dann übertragen sich die s=: 
Au Invarianten und funktionellen Abhängigkeiten. Hieraus und aus den weiter oben genannten  ...) 
E} Gründen erscheint es als am zweckmäßigsten, sowohl 'B als auch 3 kovariant-kontravariant 
zu definieren. Aus diesem Grunde ist AUCH in $3 diese Varianz bei der Definition von ® be- 

 tont worden. 
9) H. Richter, Dasisotrope Elastizitätsgesetz. Z. angew. Math. Mech,, Bd. 28 (1948), 5.205—209. ER 
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Von Rudolf Zurmühl in Darmstadt 


I iglichkeiten einer Rechnungsabkürzung zu einer wesentlichen Herabsetzung des 
aha näher den sonst üblichen Muhoden führt, Be en 
Be der aschinenrechnung 


zientendeterminante und der Kehrmatrix, und es vereinfacht sich wesentlich bei symmetrischem 
schema. i 

A melhod Jor the solution of systems of linear equations is described that is till now widely unknown and 
generally a ale, Combini He caloulus of matrices and the ae he work of recko- 


ning based on the use ofa machine one ; the computing and 
unking work, compared Zi the a ehe A reckomii dcheme wii added sum proofs and certain rules 
for Ihe practical use of the. calculating machine are given, and illustrated by a numerical example. The method 
is aypt ieh for the computation of both the determinant. of the coefficents the reciprocal matrix. It is much 
simplified, if Ihe matrix of the coefficenis is symmetrical. 

L’auteur decrit une methode pour la solution d’un systöme d’öquations lintaires jusqu’ü maintenant 
presque inconnu et d’application generale. En combinant le caleul de matrix et "usage d'une machine & cal- 
culer, on arrive & gagner beaucoup de temps. Un schema et certaines rögles sont donn£s et illustres par un 
exemple numerique. ‚La meihode peut servir pour le calcul et des determinants des coefficients ei de la matrix - 
reciproque. La matrix des coefficienis #ant symötrique, le proctde se simplifie. 

B sroü pabore ommceHBaeTca 10 CUX UOP Mao H3BecTHHÄ, Bcerna HPHTONHBIÜ METOX 
PeIIeHHA CHCTEM AMHEeÄHBIX YpaBHeHnK. Bıaropapa KOMÖHHAILUH TeOPHH MATpuI H BO3MOK- 
HOCTeÜ YMeHbBINCHNSI BEIYUCHHTENBHOH PAOOTBL IPH HCIONB3OBAHHH APHhMOoMEeTpa HTOT MeTon 
IPMBONUT, IIO CPABHCHHIO C APYTUMH MeToNaMnu, K SHAYHTEIBHOMY COKPAIeHHIO BEIYHCAHTEIB- 
HOK paboTs :u BaumcHu NPOMemyTOYHBIX pesyAısTaroB. JlpumeHureisHo u 9TOomy MeTony 
AAMTCH PaCyeTHaA CXeMA C KOHTPONEM CYMM H HEKOTOPLI® IPABHAIA IIPAKTHYECKOTO BEIJUCHEHHH 
IIPM NOMOLH ApuPMOMeTpa U BCe 9TO MWILIOCTPHPYeTcH Ha YHCHEHHEIX IPHMepax. HTOT MeTON 
HpHrOoNeH Taxe MIA BEIYUCHECHUA ONpeneuuTeneä cCHcTeMm ypasHeHnnf u OÖpaTHBIX MATpun,. 

a OH 3HayHTessHo JIIPomaeTcH B CIyyae CHMMETPHYHOH CXeMBI KOShHIHEHTOR. 


1. Einleitung 

Die mathematisch durchaus elementare, praktisch jedoch fast immer sehr mühsame Auf- 
gabe der Auflösung umfangreicher linearer Gleichungssysteme hat bekanntlich durch” Gauß 
in Gestalt ds Gaußschen Algorithmus eine Lösung gefunden, welche bis heute 
mit Recht als der vollendete Auflösungsweg linearer Gleichungssysteme angesehen worden ist, 
sofern nicht, was hier außer Betracht bleiben soll, überwiegende Hauptdiagonalglieder der 
Koeffizientenmatrix eine iterative Behandlung des Gleichungssystems nahelegen. Trotz der 
hier erreichten weitgehenden Schematisierung aller Rechenvorgänge und der durch systemati- 
schen Einbau laufender Kontrollen erzielten Sicherheit zählt die Auflösung umfangreicher 
Gleichungssysteme indessen auch heute noch zu den gefürchteten Rechenproblemen. Es hat 
daher nicht an zahlreichen zum Teil höchst geistvollen Vorschlägen gefehlt, den Arbeitsaufwand 
zu vermindern. Ein kritischer Vergleich dieser Vorschläge mit dem Gaußschen Algorithmus 
zeigt jedoch immer wieder, daß das G au ß sche Vorgehen allen anderen hinsichtlich der Anzahl 
der erforderlichen Rechenoperationen überlegen ist !). Bei der bekannten, von Gauß in seinen 
Veröffentlichungen geübten Zurückhaltung kann man vermuten, daß Gauß selbst bereits 
viele dieser Versuche gekannt, ihnen aber den von ihm angegebenen Weg vorgezogen hat. So 
ist auch schwerlich zu erwarten, daß die Anzahl der zur Auflösung erforderlichen Operationen 
noch nennenswert herabgesetzt werden kann. Ein Fortschritt in der Lösung der Aufgabe ist 
daher wohl nur in anderer Richtung zu suchen. 

Im Jahre 1938 ist von dem polnischen Astronomen Th. Banachiewicez ein Ver- 
fahren veröffentlicht worden ?), welches eine entscheidende Herabsetzung des Arbeitsaufwandes 


!) Vgl. etwa O. Heck, Über den Zeitaufwand für das Berechnen von Determinanten und für das 
Auflösen von linearen Gleichungen. Diss. T. H, Darmstadt 1946. Vgl. aber auch insbesondere den mir erst 
während der Drucklegung zur Kenntnis gelangten zusammenfassenden Bericht: E. Bodewi & ‚ Bericht übert 
die verschiedenen Methoden zur Lösung eines Systems linearer Gleichungen mit reellen Koeffizienten. Proc, 
Akad. Wetensch. Amsterdam 50 (1947), 8. 930941, 1104—1116, 1285—1295; 51 (1948), S. 53—64, 211—219. 
Die dort angeführte und gegenüber G a u ß ungünstig beurteilte Methode von Boltz-Banac hiewicz 
zur stufenweisen Berechnung einer Kehrmatrix hat mit dem hier zu beschreibenden Verfahren nichts zu tun. 

”) Th. Banachiewiez, Möthode de rösolution numörique des equations linsaires, du caleul des 


döterminants et des inverses, et de röduction des formeg quadratiques. Bull, intern. Acad. Pol i., De 
a Naben, q’ \ olon. Sci., Ser. A, 


NP oben bestimmte Regeln . B # r 
und durch ein Zahlenbeispiel erläutert, Das Verfahren eignet sich gleicherweise zur Berechnung der Kocffi- 
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22 allenfalls liegende Vorzug wird freilich durch den schwerwiegenden Nachteil erkauft, daß die 
Multiplikation nicht mehr assoziativ bleibt. 


und zwar ohne eigentlich zwingenden Grund, einer von ihm aufgestellten besonderen Matrizen- 
symbolik, der von ihm sogenannten Krakovianen oder Krakauer Matrizen %) bedient hat. Es ? 


. 


entstand so der Eindruck, als sei die Anwendung der neuen Rechenvorschrift an die Erlernung 
eines besonderen Kalküls von umstrittenen Vorzügen ®) gebunden. Die Krakovianen unter- 
scheiden sich von den gewöhnlichen‘ (C a yle y schen) Matrizen dadurch, daß ihre Multiplikation 
nicht in der üblichen Form Zeile mal Spalte, sondern als Spalte mal Spalte erfolgt. Der hierin 


Sieht man indessen von dieser für den Rechnungsgang ganz unwesentlichen Äußerlich- 
keit ab, so liegt der Kern des neuen Verfahrens 1. in der Heranziehung des Matrizenkalküls und 
2. in der Ausnutzung der in der Rechenmaschine liegenden Möglichkeiten der Rechnungs- 
beschleunigung durch das Auflaufenlassen mehrerer Teilergebnisse zum Gesamtergebnis ohne 
Ablesen und Niederschreiben der Zwischenwerte. Erst aus der Vereinigung dieser beiden moder- 
nen Hilfsmittel resultiert, trotz im wesentlichen unverminderter Anzahl der Elementaroperatio- 
nen, eine so beträchtliche Ersparnis an Arbeits- und vor allem an Schreibaufwand, daß das ' 
Verfahren wohl als ein ganz entscheidender Fortschritt in der Behandlung linearer Gleichungs- 
systeme seit dem Gaäußschen Algorithmus angesehen werden darf. Im folgenden sei das 
Verfahren in der Sprache der üblichen Matrizen wiedergegeben, wobei zugleich ein mit Summen- _ 
proben versehenes Rechenschema sowie eine Reihe von für die Zahlenrechnung nützlichen Hin- 
weisen angegeben wird. Im Falle symmetrischer Koeffizientenmatrix reduziert sich der Arbeits- 
bedarf auf etwa die Hälfte. 


2. Prinzip des Verfahrens 
Gegeben sei ein lineares Gleichungssystem 
411% FA +. FA m — 4 


yatea tt. Fam 1.0... rn. (1) 
Anı%n + ima%z H--- I Yn Cm = In 
mit nicht verschwindender Koeffizientendeterminante A = |a;, |, wofür wir in Matrizenform 
kürzer schreiben können: 
| a BEL N Fe Be FA (1) 
mit den Matrizen 
41 &ı2 --- Yın ? 2 a i 
= 1: @ı Qgs > dan f) | ” |, era en ee! (2). 
Anıtng- : - Ann In Ay, 


Führen wir noch die um die rechten Seiten a; erweiterte Matrix X = (, a) und den um die 
Komponente —1 erweiterten Vektor % ein: / 


Ar » + Ayn lı X 
EA Dr 1 LEN ea Fe Re; 
Ayı-:: Ann 4y TE 
—1 
so schreibt sich das Gleichungssystem (1) in der homogenen Form: 
Te Br ER 1 


Bei der etwa nach dem Gaußschen Algorithmus durchgeführten fortgesetzten Eli- 


-mination wird. das System (1) in,ein sogenanntes gestaffeltes System 


bu +1 +--- +bın&a=bı | 
bat, ht... td m—b; 


. A I, — b, 


3) Nachdem ich durch Zufall und ohne Autorangabe von dem Verfahren Kenntnis erhalten hatte, war 
mir Herr Brofesshr Walther, Darmstadt bei der Auffindung der Quellen in dankenswerter Weise behilflich. 
4Th.Banachiewicz, La trigonomötrie spherique et les voies nouvelles dans l’astronomie math6- 
matique.. Bull. Intern. Acad. Polon. Sei. Ser. A, 1929. 

5) Vgl. Bodewig, Zbl. Math. Grenzgeb. 17, 8.416. 
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' mit sich bringt und darum größte Beachtung verdient hätte ®). Leider ist dies keineswegs dr 
‚Fall gewesen, hauptsächlich wohl dadurch bedingt, daß sich der Verfasser bei der Darstellung, 
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Die Koeffizienlenmatrix « des Gleichungssystems (1) wird also in eine sogenannte Drei 
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eeksmatrix ® überführt, deren Elemente unterhalb der Hauptdiagonalen aus Nullen Be: 
stehen. Aus dem gestaffelten System (4) können dann die Unbekannten x; der Reihe nach be- % 
rechnet werden, anfangend mit x, aus der letzten, =,—-, aus der vorletzten. ‚Gleichung usf. | 
Der Grundgedanke des Verfahrens von Banachiewicz besteht nun darin, de 
Matrix X anzusetzen als das Produkt zweier spezieller Matrizen ® u € in der Form 


2 4 RA ee A 


bzw. jr ei | 
Y— c®» BEN a a ARTE be ER RER” (6°), . e j 
wobei B bzw. B im wesentlichen die in (5) angegebene Form der Dreiecksmatrix haben sol, 


“ während die Matrix € derart gewählt wird, daß die Elemente 5,, und c;; schrittweise einzeln Mi i 


aus der Forderung (6) ermittelt werden können. Die Matrix ® wird noch dahingehend ver- 
einfacht, daß ihre Hauptdiagonalelemente b;; = 1 gesetzt werden, wodurch sich die Aufrechnung 
der Unbekannten aus dem gestaffelten System (4) erleichtert, indem die sonst notwendigen 
Divisionen durch die d5,,; fortfallen. Die Matrix €, welche nur die Rolle einer Hilfsmatrix zur 
Bestimmung der allein interessierenden Elemente b;; (* < k) spielt, wird nun gleichfalls in Form 
einer Dreiecksmatrix angesetzt, wobe; jetzt die Elemente oberhalb der Hauptdiagonalen 
aus Nullen bestehen. Die Zerlegung (6’) lautet damit ausführlich: 


a le nd, 2020,29 RAR RAN 
(ig Aga ... Agn a 651 Ea9 0 Re 01 Bag... Dan da le (6''). 


Anı@ny... Ann An) Cry na cna ir OEL OR 


Nach einem bekannten Satz der Determinantenlehre gilt mit (6) für die Determinanten A,BC 
der Matrizen W, ®, € die Beziehung A=CB oder, da die Determinante von ® offensichtlich 
ee x 

1 En en ER RE ae ER ie (7). ii 


Man erhält also, gleichsam an Nebenprodukt, aus den Diagonalelementen der Hilfsmatrix e 


7 
den Wert der Koeffizientendeterminante A. \ 


Da die Koeffizientenmatrix X als nichtsingulär vorausgesetzt wird (A =+ 0), so läßt sich 
das Gleichungssystem (1) stets so anordnen, daß die Matrix € die gewünächte Form erhält. 
Erforderlichenfalls kann man dies durch Umordnen der Spalten von X (Umnumerieren der Ei 
Unbekannten «;) herbeiführen, vgl. hierzu Nr. 8. Einstweilen sei angenommen, daß die Matrix € | 
von selbst die Dreiecksform mit e;; <= 0 annimmt, ohne daß hierzu eine Umordnung der Aus- | 
gangsmatrix notwendig ist. | 


Daß die abgeänderte Matrix B das Gleie 'hungssystem mit den BEER Unbekannten % 1 
erfüllt, folgt aus (1) mit (6’). Danach ist 4 
AE=-CBi=0, | 

und da die Matrix € mit nichtsingulär ist (C + 0), sg, folgt hieraus: | 
BT=0 bzw. BEIN, RE 


3. Bestimmung der Elemente b;, und c;, 


Der besondere Bau der Matrizen ® und & erlaubt eine schrittweise Bestimmung der Ele- 
mente bj, und c;, aus der Forderung (6). Die Multiplikation zweier Matrizen C und ® in der 


Reihenfolge EB erfolgt bekanntlich so, daß das Element %;, der Produktmatrix W = CB gleich. 
dem skalaren Produkt deri-ten Zeile des ersten Faktors & mit der k-ten Spalte 


des zweiten Faktors ® gesetzt wird, wobei die Reihenfolge der Faktoren im allgemeinen nicht 


a Hr 72 


M hu - ie 


sprechenden Element a;; der gegebenen Matrix, vermindert um das skalare Produkt aus alle A 
oberhalb 5; » ‚befindlichen. Elementen d,, mit den entsprechenden Elementen c;, der Be 
7% treffenden Zeile , das Ganze noch a durch das Diagonalelement c;; dieser Zeile. Ein 


Er 


a ae are rn Aa 
204 


8 A 6; fürö—=k,k+1,...,n. Auflösen der Gl. (9) nach der jeweiligen a £ 


al [ 
I h Be auch der leineren, der beiden Zahlen ü ; oder k bzw. 2: 
= fi “ e Pr % r 
rechne nun ‚(mit Rücksicht aut eine er erforderliche Spaltenumordnun 
gangsmatrix) zweckmäßig spaltenweise, beginnend mit der ersten Spalte von Y und. 

durchläuft jede Spalte von oben nach unten. Dabei tritt dann in jeder der En aa 
genden Gln. (9) jeweils ein neues, also unbekanntes Element: Dir oder Ci, auf, und zwar in 


der k-ten Spalte von ® und & zunächst die Elemente. bi, füri—=1, 2, ...,6—1, danach die 


bin an ü u & oder Cr wur >> k) ergibt die er % RR 


ea Be; 
z a dene >= k') et (10b). 
De Vorschrift lautet: Ein-Element 5;, auf dem Platz ;, ki < k) berechnet sich aus dere Say i 


me ‘1 dh E Er 


= (4, bı Kr fie br. 


= Sad) »— ip ba... 


' Element c;, auf dem Platz i,k (6 >%k) berechnet sich analog, aus dem entsprechenden alten 2 
Element a;,, vermindert um das skalare Produkt aus allen vor cy; befindlichen Elementen 
€;, mit den entsprechenden Elementen box der betreffenden Spalte k. Wegen b;;—= 1 erübrigt 
sich hier eine Division. Die Vorschrift (10) gilt auch noch für die Spalte der reehten Seiten @; 
und 5;, wenn man formal setzt a; — @;, n+1 5 = bi,nzı. Sie gilt auch noch für eine unter 
X und € anzubringende Kontrollzwecken dienende Summenzeile —s; und —o,, wenn man 
setzt — 5, = An4ı,% —0;: = Cn+ı,m Vgl. den folgenden Abschnitt. 
Für die erste Spalte der &-Matrix folgt sogleich: 


Berge ne 


Für die erste Zeile der 3-Matrix gilt entsprechend: 
Par=aziau|' (1,2, :..,n +1) TRETEN FE 


4. Rechenschema und Proben 
Für die praktische Durchführung der Zahlenrechnung vereinigen wir die beiden Dreiecks- 
matrizen ® und € zu einem einzigen rechteckigen Zahlenschema, wobei in der Hauptdiagonalen h 
nur die Elemente c;; eingetragen werden, während die nicht weiter benötigten Elemente b;;— 1 Be, 
fortgelassen werden können. Dieses Schema, das wir fortan kurz die „,BC-Matrix“ nennen +, 
wollen, schreiben wir etwa unter das Schema der gegebenen Koeffizientenmatrix Aund erhalten 


„Tadel 1. Rechenschema zur Auflösunglinearer Gleichungssysteme | Er = h 
nach dem, Matrizenverfahren 


N Caı Eng Eng he # Can by | En 

— 0% | 00 en | —=0 | lg 
Probe ? ? ? ? = Zaun 
k=1: 2 3 n | b; vi 
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denkt man sich die Matrix X einfach um eine weitere Zeile mit den s, erweitert, so erhält man 


auch in der Hilfsmatrix C eine weitere Zeile mit den entsprechenden Elementen o,. Diese Größen 
müssen dann, wie man sich leicht überlegt, gleich den Spaltensummen der &-Matrix sein, also 


0% == C% +. rt: + nk Fe A De er 3 . (13b). 
Praktisch verfährt man besser, ähnlich wie bei Gauß, so, daß man die negativen Werte —s; 


. und —o;, benutzt und zur Kontrolle die Gesamtsummen 


% Gr FH, rt. mat)... 2.2.20. 
bildet, welche Null ergeben müssen. Infolge unvermeidlicher Abrundungsfehler können davon 
abweichend wenige Einheiten der letzten Dezimale auftreten, die man, in Einheiten dieser Dezi- 
malen, in einer letzten Kontrollzeile laufend aufschreibt. Im Rechenschema der Tafel 1 ist zu 
beachten, daß die Spaltensummierung dabei nur über die C-Matrix, also die Glieder auf und 
unter der Hauptdiagonalen zu erstrecken ist. Dem gleichen Kontrollvorgang werden übrigens 
auch noch die Spalten der rechten Seiten unterworfen. Das der a;-Summe s entsprechende 
Element o muß dabei Null ergeben, da hier keine c--Elemente mehr auftreten. In Tafel 1 sind 
die Kontrollplätze durch Fragezeichen angedeutet. 

Sind sämtliche Elemente der B@-Matrix bestimmt, so beginnt de Aufrechnung 
der Unbekäannten x; aus dem gestaffelten 'Gleichungssystem (4) mit 5;; = 1. Man er- 
hält, anfangend mit der letzten Unbekannten «,, der Reihe nach aufsteigend: 


En Fe re )° 


2x, = b; — In Din Un b,, ee T+ b;, i+1 


Sind die n—i letzten Unbekannten bereits bestimmt, so errechnet sich die i-te x; demnach aus 
der rechten Seite b;, vermindert um das skalare Produkt aus allen schon vorhandenen Unbe- 
kannten mit den Elementen d;; der i-ten Zeile. Man durchläuft die Spalte der x, und die Zeile 
der b;, gleichsam in umgekehrter Richtung wie sonst und hört dabei wieder, wie stets, vor dem 
Diagonalglied der 8C-Matrix auf. 

Sind schließlich alle Unbekannten berechnet, so folgt noch die Probe durch Ein- 
setzen der rn Unbekannten &,; in das Ausgangsgleichungssystem gemäß 


Yıkı Füst3 t... + m; = 0 En a ES 3 2 LEE (16), 


wobei? = n -+1sich auf die letzte Summenzeile der A-Matrix bezieht. Man bildet also wiederum 
skalare Produkte aller Zeilen der A-Matrix mit der um die Komponente —1 verlängerten Er- 
gebnisspalte der &;. Die Ergebnisse der Proben, also Nullen bzw. wenige Einheiten der letzten 
Dezimalen, trägt man in die der A-Matrix angefügten Probespalte ein. Bei ungünstiger Koeffi- 
zientenmatrix (relativ sehr kleine Koeffizientendeterminante) können die Einsetzungsreste 
übrigens gelegentlich auch beträchtliche Werte annehmen, ohne daß die Rechnung falsch sein 
müßte; doch ist dies ein seltener Ausnahmefall. 


5. Die Zahlenreehnung 

Zählt man die beim Matrizenverfahren erforderlichen Einzeloperationen, also die Anzahl 
der Additionen, der Multiplikationen und Divisionen ab, so kommt man ungefähr auf die gleiche 
Anzahl wie beim GauBschen Algorithmus. Wollte man also etwa die Multiplikationen einzeln 
vornehmen, die Ergebnisse aufschreiben und alles addieren, so hätte das Verfahren gegenüber 
dem bewährten Gaußschen Vorgehen nicht den geringsten Vorteil. Eine Arbeitsersparnis 
ergibt sich erst dann, wenn man die für das Verfahren charakteristischen skalaren Produkte 
unmittelbar mit einer Rechenmaschine, d.h. durch Auflaufen- 


*) Die Summenproben finden sich in der unter 2) zitierten Hauptarbeit noch nicht, wohl aber in an- 
derem Zusammenhang in einer weiteren Arbeit: Th. Banachiewiez, La rögle de Chid, cracoviens et 
matrices. Bull. intern. Acad. Polon. Sci. Ser. A, 1938, S. 405-412. 
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RE inze rn p rodu kt e. Nurin EN Verbindung des Fer an es a durch die Eig 
‚art der, Rechenmaschine bedingten Rechnungsabkürzung kann die entscheidende u: 
heit des Matrizenverfahrens ‚gegenüber den bisherigen Rechenmethoden gesehen werden. 

"Verfahren steht und fällt mit. der sachgemäßen Benutzung. ‚der Rechenmaschine.. Demg ei 
‚ist die Beachtung einiger einfacher Regeln und Kunstgriffe für ein rasches und sicheres Arbei 
wesentlich, die indessen leicht zu en sind, so daß auch der im Maschinenrechnen Ungeü 
sich bald hinreichend einarbeitet. ' 
° Zunächst ist es mit Rücksicht of &itheitliche Kommastellung Asa, die Fichlen ER 
| SER falls notwendig, durch Multiplizieren mit einer geeigneten Zehnerpotenz auf die gleich 
Br Größenordnung wie die übrigen Koeffizienten Q;, zu bringen gemäß a; = 10”a;. Das bedeutet. 
den Ersatz der Unbekannten x; durch neue Werte 2; = 10”, welche von der "Größenordnung En; 


der I sind. Alle Elemente der X-Matrix erhalten dann die gleiche Stellenzahl nach dem Komma. 
Wie die GIn.(10) und insbesondere (11) und (12) erkennen lassen, sind die Elemente Ei 
von der Größenordnung der Koeffizienten a;,, die Elemente bir dagegen von der Größenordnung. 
der 1. Die beiden Zahlenarten der ®C-Matrix werden also im allgemeinen nicht mit gleicher 
5 Stellenzahl nach dem Komma, sondern, mit Rücksicht auf -volle Ausnutzung der Genauigkeit, 
F er gleicher Anzahl geltender Stell en. BUOChBEn 2 B. durchweg mit 6 höcr 7, ja rr 
gl. das Zahlenbeispiel in Nr. 6. ei 
Mit Rücksicht auf die bei den bi; erforderliche Division empfiehlt es sich, bei a Ai Bu: 
e; rechnung der skalaren Produkte in (10a) und (10b) die Faktoren c;, grundsätzlich im Ein- 
 „stellwerk (Tastatur, Einstellhebel) der Maschine einzustellen, die Faktoren d,, hingegen 
dann wie üblich indas Umdrehungszählwerk laufen zu lassen. Das erste Glied a;; ER 
.- » wird wie ein c-Element behandelt, multipliziert mit dem Faktor 1. Die Ergebnisse e;, werden 
\ ' dann unmittelbar im Ergebniswerk der Maschine (im Maschinenwagen) abgelesen. Für die je 
Ergebnisse b;, ist der gleichfalls im, Ergebniswerk erscheinende Klammerwert der Gl. (10a) 
noch durch das Diagonalglied c;; zu dividieren, was durch Eintasten des c;;-Wertes wie bei den 
anderen. c-Werten und anschließende Ausführung der Division geschieht, Das Ergebnis b;, 
erscheint dann in derüblichen Weise im Umdrehungszählwerk der Maschine. Bei Beachtung dieser 
Regeln arbeitet man mit fester Kommastellung bei der gesamten Berechnung der B&-Matrix. 
Die Produkte werden mit den in den Formeln (10) angegebenen unter Beachtung der den 
‘Zahlen eigenen Vorzeichen gerechnet, wobei sich negative Produkte selbsttätig abziehen. Nega- 
tive Ergebnisse bei den c;, erscheinen dann als Komplementzahlen, z.B. ...9943,8752 = 
—56,1248,. Wird der Klammerwert bei b;; negativ, so wird die Division nicht wie üblich durch 
Abziehen, sondern durch Hinzulaufenlassen des Divisors’c;; bis auf. die Null aus- 
geführt bei positiv (nicht wie sonst_bei der Division negativ) zählendem Umdrehungs- 
zählwerk, wobei freilich eine automatische Division nicht mehr, wohl'aber noch ‚‚Stoppdivision‘ 
möglich ist. Das Vorzeichen von b;; richtet sich dann natürlich noch nach dem von c;,. 

‚ Die Ergebnisse &; sind von der gleichen Größenordnung wie die b;,. Bei der Bildung der 
skalaren Produkte in (15) müssen also die Kommaeinstellungen auf der Maschine u. U. abgeändert nt 
werden. Die Probe nach (16) erfolgt dann wieder mit der alten Stellenmarkierung, wobei die 2; 

4;, in das Einstellwerk, die &; in das Umdrehungszählwerk gegeben werden. ; 
An Rechenmaschinen eignet sich grundsätzlich jede der üblichen nach dem Additions- 
prinzip arbeitenden Maschinen. Doch ist bei Maschinen, welche nicht über eine durchgehende 
 Zehnerübertragung im Ergebniswerk verfügen (Staffelwalzensystem wie Rheinmetall-Borsig 
oder Archimedes), Vorsicht geboten, da hier das Ergebnis beim Auflaufen der skalaren Produkte 
durch die auf der Vorderseite des Maschinenwagens aussetzende Z.ehnerübertragung leicht ver- 
ER fälscht werden kann, was indessen bei einiger Aufmerksamkeit im gegebenen Augenblick be- 
’ merkt wird. Alle übrigen Maschinenbauarten (Sprossenradsystem wie Brunswiga, Zahnstangen- 
‚system wie Mercedes-Euklid, Schaltklinkensystem wie Hamann) sind ohne weiteres und mit 
völliger Sicherheit verwendbar. 


6. Ein Zahlenbeispiel 
Gegeben sei das lineare Gleichungssystem für 4 Unbekannte: 

43,570, — 23,84, — 51,653 - 19,44, = 0,7362 

62,89%; + 84,97, + 21,842, — 39,352, = 1,1872 

37,48x, + 93,24; + 84,395 + 26,752, = 0,3875 

-19,37x, + 54,382, + 14,598, + 62,85, = 0,5738 
Die rechten Seiten werden durch Multiplikation mit 100 auf die gleiche Größenordnung wie die 
Koeffizienten gebracht, man reehnet also mit abgeänderten Unbekannten x; = 100x,. Nach 
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NEN Schutzstellen Se die ROChBEIn ee sechs- b 
In Zahlentafel 1 ist die gesamte Zahlenrechnung enthalten. Die Zahl fr. it 
| und dritten Spalte entsteht beispielsweise aus f 


Sa — 84,3900 — 37, 4800 ae | ‚185449 — 113,7477- : 0,807438 = or. . 
ei Die en Ps in der dritten Zeile der d;-Spalte entsteht aus % F 
Be . b,= (38,7500 — 37,4800 : 1,689695 — 113,7477 - 0,104330) : 36,9764 = = _0,985685. 


7. . Zahlentafel 1. "Zahlenbeispiel zur Auflösung linearer Gleichungssysteme nach dem 
Bu.’ Matrizen verfahren 
3 : kl 2 3 4 fi Probe + 10% 


1 43,5700  —-23,8400  —51,6500 
= 2|  62,8900 _ 84,9700 21,8400 
Ei’ DE RR 37,4800 93,2400 84,3900 
N 4| —-19,3700 54,3800 14,5900 


— sr | —124,5700 —208,7500 —69,1700 


N ve 43,5700 | 0,547165 —1,185449 SR as 0,637671 
| ’ 2| 628900 | 119,3812 0,807438 —0,564663 | 0,104330 | 1,404457 

2, \ | 3| 374800  113,7477 | 36,9764 2,008209 | —0,985685 || —1,448886 
we: \ 4| —19,970 43,781 4,7229 | 184,0189 0,230654 | 0,230654 
E — 0% F-124,5700 —276,9003 [6,7465 —184,0189 0,0000 | —ı 
% Probe 0 0 0 0 j 
| | ra! 2 3 ur B; 2; | 
P Die Ergebnisse sind, auf vier Stellen = nen entsprechend den vierstellig gegebenen Aus- 
3 “ gangswerten: x 
{ | | 2 =  0,006377 

=  0,014045 

& = —0,014489 ° 

<=  0,002307 


7. Vereinfachungen bei symmetrischer Matrix 
Wie bei den meisten übrigen Verfahren zur Gleichungsauflösung ergeben sich auch beim 
Matrizenverfahren Vereinfachungen im Falle symmetrischer Koeffizientenmatrix ?). Dabei ist 7 
nämlich: 


a . 


d.h. auch die 8C-Matrix ist in gewisser Hinsicht symmetrisch, wenn man auf der B-Seite, also 
oberhalb der Hauptdiagonalen noch die Division durch die Diagonalelemente c;; berücksichtigt. 
Bei der Berechnung eines Elementes c;, wird demnach das gespiegelte Element b;; der BC- 
Matrix einfach dadurch mitgeliefert, daß man das Ergebnis c;; noch durch das in der gleichen 
Spalte darüber befindliche Diagonalelement c;; dividiert. Eine besondere Berechnung der b- 
Elemente erübrigt sich somit mit Ausnahme der rechten Seiten b;. 

Die Richtigkeit dieser Behauptung ergibt sich durch Induktionsschluß wie folgt unter i 
EESHELSNE der Symmetriebedingung a;, = a;;. Für die erste Spalte und Zeile der 8€-Matrix 
gilt: 

G1>%ı 
di tu; >=, =: 


Für die zweite Spalte und Zeile liefern die Gl. (10) hiermit:, 


se 
a B 21 
a= ta — bad —i,— 
1 
Cds dr Cd Be 
22 ai; = Ai 21 1877 ET Ne 
11 


‘) Hierauf hat Herr Dipl.-Math. Gerh.Schreier, Darmstadt aufmerksam gemacht. Bei Ba- 
nachiewiez selbst findet sich darüber nichts. 
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‚zeichen umkehrt, so berechnet sich die Koeffizientendeterminante jetzt zu 


&s 
iz = Ca ust. 


MN 3: 
hemmt en bar = a, ‘= EN: 
23) 
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Ce ER EZ 8. Abweichungen v vom Regelfall hier 

Das Verfahren. in der bisherigen Form ohne De versagt, sobald eines & 
> Disgonalelemente c4 = 0 wird, was genau dann eintritt, wenn eine Hauptabschnittsdeterminant 
‚der Ausgangsmatrix N verschwindet. Da andrerseits A= C =# 0 sein soll, kann der Fallc;,—=0 \ 
nur in der Weise eintreten, daß die Spalten von & gegenüber der normal angenommenen Drei- 
ecksform vertauscht sind. Das aber kann dann durch ein entsprechendes Vertauschen in der 
Ausgangsmatrix A behoben werden. Im Rechenschema, Tafel 1, bedeutet dies lediglich ein 
Abändern derNummern k am Kopf des Schemas, wobei jetzt freilich die abgeänderten Nummern 
mit den entsprechenden unveränderten Nummern k am Fuß des Schemas in Verbindung zu ‘ 
bringen sind. Sobald also im Laufe der Rechnung in einer k-ten Spalte das Diagonalelement 
C+ = 0 wird, vertauscht man die beiden Nummern k und k --1 am Kopf der X-Matrix und 
“rechnet jetzt die k-te Spalte der 8C-Matrix mit der (® + 1)-ten (jetzt mit k bezeichneten) Spalte 
_ der X-Matrix, und ebenso die (k -+ 1)-te Spalte a 8C-Matrix mit der k-ten (jetzt mit k +1 
bezeichneten Spalte der X-Matrix. Dabei müssen natürlich die bereits berechneten b;,-Werte 
oberhalb des Elementes c,, ausradiert und neu gerechnet werden. Die Umnumerierung muß 
dann zum Schluß auch bei der Reihenfolge der Unbekannten berücksichtigt werden. Ver 
schwindet nach der Spaltenvertauschung das Diagonalelement noch einmal, so ist die Vertau- 
. schung zu wiederholen, d.h. die k-te Spalte der X-Matrix ist durch die (k + 2)-te Spalte zu er 
setzen usf. Da bei einer Vertauschung zweier Spalten die Determinante einer Matrix ihr Vor- 

ER 


1,3 


‚A=(-1C 
wenn z die Anzahl der Spaltenvertauschungen ‚angibt. 


Die Spaltenvertauschung wird man übrigens auch dann vornehmen, wenn ein RL 
element c,,im Vergleich zu den übrigen c-Werten auffallend klein wird. Auf diese Weise schützt 
man sich vor u. U. beträchtlichen Genauigkeitsverlusten. 


= (18), 


” 


Beispiel: 22, +32, +42,= 13 
4x, +62, — 220, = —4 
— 43% +2,3= 1 eh 
| 1 32 a; Probe E 
1 2 3 4 13 € Er 
2 "4 6 —2 —4 % 
3 —l 3 2 al . 
EB alu] 0 
’ a N RE 
2 u & 
3 . 
a 
f Probe 
Ergebnis: wie me. ; 


A=C-I=+%. 


9. Berechnung der Kehrmatrix 
Die Berechnung der Kehrmatrix (der inversen Matrix) A! — (&;,) zu einer gegebenen 
Koeffizientenmatrix WA = (a;,), definiert durch 
AAI=-ATA=-E (19) 


6% 
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mit der Einheitsmatrix y 


- vollzieht sich beim Matrizenverfahren in der gleichen Weise wie beim G au ß schen Algorithn 
einfach dadurch, daß man die Auflösung eines Gleichungssystems mit der Matrix A und r 
den n-Spalten der Einheitsmatrix an Stelle der einen Spalte a; als rechten Seiten durchführt, | 

Im Rechenschema der Tafel 1 treten also an die Stelle der einen e a; die n-Spalten dr 

Einheitsmatrix. Entsprechend ergeben sich an Stelle der einen Spalte b; der ®C-Matrix Das; (7 
n-Spalten. Aus ihnen berechnen sich dann n Ergebnisspalten mit n?-Zahlenwerten z;,, deren # 

Schema unmittelbar die gesuchte Kehrmatrix darstellt gemäß RE 2a R 


h L = %n »»IN 0 re r Yu a Be A; Feskäck 


{ 10. Zum Arbeitsbedarf vr AN 
9%. Willman den Arbeitsbedarf des Matrizenverfahrens vergleichen etwamitdemdes Gauß- 
schen Algorithmus, so kann man zunächst einfach die Anzahl der niederzuschreibenden bzw. 
der der Rechenmaschine zu entnehmenden Zahlenwerte abzählen. Macht man dabei auch im 
Gaußschen Algorithmus, soweit möglich, von der Abkürzung durch Auflaufenlassen yon 
Teilprodukten im Ergebniswerk der Rechenmaschine Gebrauch, so kommt man unter Einrech- 
nung aller Summenproben sowie der Aufrechnung der Unbekannten und der Einsetzungsproben 
zu folgenden Anzahlen z. Einschließlich der gegebenen Koeffizienten a;, und der rechten Seiten @; » 
sind niederzuschreiben: / 


r er beim Matrizenverfahren: zu =2n®-+7n+3=(2nr-+1)(n +3) 


- 


x 


beim Gauß- Algorithmus: 2a (2m 415n +19) — r (zu +8n +16). 


die Anzahlen 2’ der der Maschine zu entnehmenden Zahlenwerte. In Zahlentafel 2 sind diese 
Anzahlen für einige Werte der Zahl n der Unbekannten zusammengestellt Die Anzahlen ver- 
Be .® ringern sich natürlich bei beiden Verfahren, wenn die Koeffizientenmatrix, was praktisch häufig 
= der Fall ist, teilweise aus Nullen besteht. 
Zahlentafel2. Vergleich der Anzahlen der  niederzuschreibenden 


bzw. der der Rechenmaschine zu entnehmenden Zahlenwerte beim 
Matrizenverfahren (z,,) und beim Gaußschen Algorithmus (24) 


e Br. Zieht man beide Male die Anzahl n (n +1) der gegebenenWerle a;, und a; ab, so erhält man | 
“ 


x 

Zahlder Niederzuschreibende Der Rechenmaschine zu : N 

Unbekannten / Zahlen entnehmende Ergebnisse j 

' 
3 42 41 30 29 

4 63 74 42 54 R 

5 s8 ‚490 58 90 
10 273 615 163 505 

20 943 3 730 523 3 310 

50 5.353 48.075 2 803 45 525 } 

100 20 7083 358 650 10 603 348 550 3 

x 


Das Verhältnis des tatsächlichen Arbeitsbedarfes ist natürlich keineswegs so kraß, wie es 
nach dieser eindrucksvollen Zusammenstellung den Anschein haben könnte, da die Berechnung 


jedes Zahlenwertes beim Matrizenverfahren im Mittel eine wesentlich größere Anzahl von Einzel- 1 
operationen verlangt als bei Gauß. Verläßliche Angaben hierüber sind indessen nicht ohne 
eine sorgfältige Analyse der Elementarzeiten möglich, unterstützt durch umfangreiche prak- 


tische Rechenversuche etwa im Sinne der unter 4) zitierten Arbeit von Heck. 

Ein nicht zu unterschätzender Vorzug des Matrizenverfahrens besteht übrigens noch in 
der völligen Gleichförmigkeit aller Rechenvorgänge, nämlich der fortgesetzten und ausschließ- 
lichen. Bildung skalarer Produkte, wobei sich auch die zur Berechnung der b;;, erforderlichen 
Divisionen durch die Diagonalelemente gleichsam als letztes Glied dieser Produktbildung glatt 
in den übrigen Rechnungsablauf einfügen, sofern man die oben angeführten Hinweise zur Ma- 
schinenrechnung beachtet. Hieraus ergeben sich möglicherweise Vorteile für eine rein maschi- ie 
nelle Auflösung linearer Gleichungssysteme durch Rechenautomaten (Hollerith u. dgl.) 
Eingegangen: 15. 5. 48. 
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I HERR Er . Obere d tag Schranken ni eine Nullstelle 2* ei einer ‚Funktion ; u" = 
Ro we Ableitu 


Fort hr Ra 


2} LEIDEN Seine Kohvergenz. beschleuni en läßt und daß es "such bei Ni Re en ‚der, 
B, Er i an ga A Schranken a ar eh BR! 2 a0 Fl 
& | 


KATHL) re polynämes du n-idme ordre en ©, qui enferment fe) dans le. voisinage de x*, a condition ‚que la 


7 NOKABBIBACTCH, KAK MOIRET es YCKOpeHa, CXORHMOCTB JTOrO MeToNa U Kak IIyTeM HEBHAUH- | 


AT, Einleitung. Die Funktion ya Ei sei in einem Intervall I ER 2< 2 Knie ig a 


sei in J f®D (2) monoton und f(x) (e=1,2,..n) dort stetig, weiter sei n — 4423 


 fle) kann Kuh _ Mo | 
Ya (x) in der Umgebung von a* vorausgeselzi — durch die Nullstellen gewisser, fa) dort, 


schinßendennt Poly mome N-ten Grades in & erhalten. Für dieses. ‚Verfahren. wird u. a. nachgewiesen, 


X 


ar of Ihe N-th order i in ® tat are riküdig f(«) in et BR u, Knenlied, (hat the derivalive en 
is monotonuos. Among other things, it is proved howto accelerate the convergence and that ihe method 7 
3 en of x* after a slight ‚modification, even if the supposilion of Ihe monotony is not realized. 


‚On peut arriver aux limites inferieures et superieures d’un zero x* d’une ‚fonetion y= fe) ü N 


derivee Tal I (x) soit monotone. On montre comment on peut accelerer la convergence ei on prouve que Tas 
methode fournit (apıes une simple modification) des limites pour u”, meme sans ‚supposer une monntonie. 


Er Bepxuss MH HWKHAA TPauuıst Kopnst 2* pyaeuun y=f(a) .MOTYT ÖBITB. HONyyeHst B 
CIyyae MOHOTOHHOCTH 1IPOUSBOAHOH FIN (x) B orpecTHocTn w* NpH IIoMorMm KOPHeü HeKo- 
TOPEIX IOIHHOMOB CTeICHU N, BARINOYAMIUHX B HTOoH Oßutacrnu hyakummo f(x). Kpome Toro 


TeIBHOIO EIO H3MEHEHUA MOTYT ÖBITB HOIyYeHEI Kessuhee a «* Tarıke u B cayuae, a 
yCHOBHe MOHOTOHHOCTH HE BBINONHeHO.  '  , 


und besitze dort genau eine Nullstelle «*. Es seif (x) *f (25) < 0. Obere und untere Schranken 2 
für &* sind durch die Nullstellen gewisser Interpolationspolynome Yan (2) gegeben, Dabei 


% "65 >1 und y.= 9a, ( x) dasjenige (eindeutig bestimmte) Polynom höchstens ( (n—1)ten f 
Grades in x, das mit f(x) in«=z, (u =1;2) im Funktionswert und in den aufeinander- 
folgenden Ableitungen nach x bis zur (%. — I)ten übereinstimmt. (f® (z,) bedeute f(x). 
Zu diesem, in einem früheren Aufsatze!) betrachteten Verfahren der Wurzelbestimmung durch x 
Einschließung, &iner Verallgemeinerung des kombinierten Verfahrens Newton-Regula falsi, seien 


hier einige Bemerkungen ae \ 
Für den einfachsten Fall» =3 (bei a,=2, 0, =]1 bzw.0,=1,a, 2) zeigt die Darsteilüng 2) 
1 st h; 
Fa) — P2,1(®) Fee. (0 — 2)" (20); 


fix) — Pı.(8) = 3] IM T. zo) 


NE < N 
daß für alle x in J. 


entweder stets P12(2) <f(x) <Pp51(a) gilt, falls dort ständig f(x) >0 | ) 
oder u ESS SDEM nr nn fm 3 Da 


Für in J Be x) schließen also die Parabeln 2.0. p,,.(2) und 95,,(8) dort f(x) ein. 
Bei f (z,) : f(#,) < 0 sind somit durch die im Intervallinnern gelegene Nullstelle von y=P3;,1(%) 
bzw. die von y = 91, (x), die wir mit x%ı bzw. x%,g bezeichnen wollen, obere und untete 
Schranken für x* gegeben. 

Der Unterschied dieser Schranken ist (bei f’(2*) = 0) klein in h (2 h=%— ©, > 0) von 
dritter Ordnung. Durch einfache Umbildung von f(x) läßt sich, wie in Abschnitt 2 erläutert 
wird, erreichen, daß die so sich ergebende neue Abweichung ”,, 1—af, o mit der fünften Potenz 
von h gegen Null geht. Daß auch bei Nichterfülltsein der eben erwähnten Monotoniebedingung 
für f” (x): das Einschließungsverfahren sich durch Übergang zu einem geeigneten Gleichungs- 


system a durchführen läßt, wird in Abschnitt 3 gezeigt. In Abschnitt 4 geben wir Beispiele für 


die Gewinnung von Schranken für &* bei Einschließung von y = f(x) durch in y quadratische 
Polynome. Schließlich werden in Abschnitt 5 bei Zugrundelegung von FR Interpolationsstellen 
weitere Ne re hergeleitet. — Die nachstehenden Überlegungen lassen sich auch für 


n—=4,5,... durchführen. Die dabei sich ergebenden Formeln sind jedoch im allgemeinen für 


2) E. Pflanz,.Z2. angew. Math. Mech. 28 (1948), 8. 114—122. f 
°) A. Marko fi 2 Jifferetrepesahnung (deutsche Übersetzung von Th. Friesendorff und 
E.Prümm). Leipzig 1896, 


' ei 1978 3% nn. mr N . „ N BR ji E 
86 Plans, Zur Bestimmufgmellr Nulltelle 
> den praktischen Gebrauch nicht genügend einfach. Wir gehen daher ni 

und beschränken uns auf diesen Hinweis. KePR b ‘ 


ET N: 
2 r ü er 


A 


7 


Be 2. Konvergenzverbesserung. Es sei (2) :f(2,) <0, f” (2) monoton E 
dort stetig. &* sei die einzige Nullstelle von f(2) im Innern von J. Dann gilt®), bei 

R: in F und für kleine = (0,2) | | 1 2 Er 

Bi: } 5 Fu Iris Bern, | 
Br | h3 &e\? &,2 Er 2 
BE an (el 


| 1 5% | Ha, | 
mit &,ı = 8,1 — Umy .. +» Im = 5 (2, + 22), fm =f' (2m); - - - , wobei die in (3) weggelassenen BE 


RS ME Glieder von höherer als dritter Ordnung in % sind. Nach (3) wird das z*einschließende Intervall FE 
u .. «$ı— af und damit der Fehler der Näherungen für z* (für kleine h) besonders klein sich 
B', ergeben, wenn f”’ (2) = fm = 0 ist; denn dann verschwindet die Abweichung 21-2 von, 

Ku höherer als dritter Ordnung. Der Ansatz \ f Ru 

E: SHE HF), Er) ner 1): Vor 


gibt nun bei passender Festlegung der Konstanten b,, b, eine Funktion v(z), für weiche 


vo =v" (2m) = 0 ist und die wie f(x) für »=x* verschwindet®). Neben der Forderung 

Er N I 5). 

muß, um der erwünschten Monotonie von v”’ (x) in J nicht zu widersprechen, auch 

3 : \ 

; ME... Vor er (6) 

i / sein, denn sonst ändert v’’ (x) in der Umgebung von & = x. sein Vorzeichen. 

Beni . 5) und (6) ergeben für die Konstanten b,, b, des Ansatzes (4) j 
EL = — Um fm 6, fm —3Imfn - N an 


IE NEN DT A a DIES AT N A i 


Weiter erhält man [auf gleiche Wejse wie die Entwicklung (3)] für das zu v(z) nach (4) und (7) 
gehörige, neue Schrankenintervall 


° 


ee =& —& De tor = 1,2 : BT: ER a 0 8 
2,1 12 SSL 12 = 55 n ++: wert 9778: 


Um 


Wegen der Gl. (6) bis (8) sei in J /(x) fünfmal stetig nach «& differenzierbar voraus- 


‚gesetzt. 
- : 

Bei dem Beispiel [()=e— 2 at so, =0, m =L it > Kan b, = = . Damit ergibt . 
sich v. (2) = 15 e’-& >20, also v’’(x) monoton für — -; SU Te 7° 0<z<1und es ist dort 8 
1 +biE —+b,£?=#0. Nach (8) ist die Abweichung 3,1 — 27,2 = — 1,5 107% zu erwarten. 
Die Auflösung der quadratischen Gleichungen p%,(2) =0 bzw. 9,,.(2) =0 führt auf*die 
Ungleichung 3,1 = 0,69307 - - < x* = In 2 = 0,69314 - - < 0,69321 -- - = z% 3, gegenüber dem 
Ergebnis°) 0,673 < 2* < 0,703 bei der gewöhnlichen Einschließung von f(x) (also bei b, =b,—0) 
durch Polynome 9,,,(%), P3,1(8): 

Bei dem weiteren Beispiel f(x) = 2? —3(1< .x*< 2) erhält man gleicherweise durch k 


Übergang zu v(#), mit y =1, 2%, =2 1,4415... < ar < 1,4430. .., und durch Wiederholung 

des Verfahrens mit x, = 1,44, x, = 1,45 unter Beibehaltung der Werte b,, b, des ersten x 

Schrittes 1,442 249 570,9... < @* < 1,442 249 57073. .. Mit b, = b; = 0 ergibt sich bei &, = 1, E 

©4=2 hier nur, 1,425., <a*< 1,463... .. 
Zur Berechnung von b,, b, nach (7) benötigt man die Ableitungen von /(«) bis zur vierten | 


Ordnung in der Intervallmitte. Man kann jedoch anstatt b,, b, auch Ausdrücke b,, b, ver-. 
wenden, die nur die ersten und zweiten Ableitungen von f(x), und zwar an den Intervallenden, 


enthalten. Die Abweichung En — der .Polynomnullstellen, welche zu der mit bs b, ge- 


®) Siehe a.a.0.t), Gl. (24a). 

‚...‘) Um die Konvergenz bei der Anwendung des gewöhnlichen Verfahrens von Newton und der er- 
weiterten Newtonschen Formeln zu beschleunigen, wurde von Fr. A. Willers (Z. angew. Math. Mech. 
16 (1936), 5. 315—316, und 18 (1938), 8. 197—200) eine Methode angegeben, mit der unser oben erläutertes Ver- 
fahren den Ausgangspunkt gemeinsam hat. Zur Verbesserung der Konvergenz sucht man in der{ Reihen- 
entwicklung für die v-te Abweichung d, (einer »-ten Näherung für x*) nach Potenzen von dy-ı den Faktor 
3 Be, es, en zu machen bzw. zum Verschwinden zu bringen, Bei Wil- 

ers geschieht dies durch Transformation der Variablen x, bei un i i 
ee serer Methode durch Umbildung der Funktion. 


re meta, 


ne 4, er 
b «30 \ 


“2 = Be a Een a, Ki en Em . DE 
= Die Auflösung von 9) ‚nach I 2 ergibt mit den Abldreingeh I I RR “ per a RR 
ee hei BeRER Le 
die Entwicklungen k\ R Ye RL BR re > 
a ö f' . 1 1 ER ” R A Be ? ar cH, | 27 i | ) A H A ara N us. - Ki 
En: Leite et SP Es s Re hats u % 
Be xl Ina" +. 3 BO=-— +, 3 ho m iR: 5, j Be 


in (10) sind 81... gewisse, hier näher nicht interessierende lineare er üben Be 
. die in (9) auftretenden fünften Ableitungen von f. Setzt man fm, ... aus (10), unter Vernach- 
BASIEUNG der S,, in m ein, so erhält man die gesuchten Ausdrücke GE a 


EN TIERE ST DER 
AREA z h [2786'2— 25 (96 o.40) 6") + (802 —+409”2)]' 
Bi: % ee 2762 —36hö’ 0" + h2 (90218 8”2)] BER 
EN Be GM [27 62 — 2h 490 0" —+40’0") + A? Ge?+46'%] 
Für ee f(x) vom Grade <A ist (wegen Sa el) Denn Den 
Der Zusammenhang zwischen b, und b,, b, und b, läßt sich für kleine A darstellen durch 
b,=b, Rp... eb, + pH... 


wobei die (für kleine Ah beschränkten) Ausdrücke 9;, y, die Funktionswerte-fi, f, far J- und 
die 8,, ... enthalten. 


. 
Für die Nullstellen z3;,, eis 0 bzw. &,, E; 2 der zu v(z) gehörigen Einschließungspolynome \ 
Pa» Pı,2 folgt dann- 


Ei a ei » En Fi pl h-(&=)7 Hr (B% +0) (12): 


TER UTEEE WEBER) 


5! Um 


Durch Heranziehung der Funktionen v (x), d (x) anstatt f(x) und bei Einschließung von v (x), v(x) ;s. 
durch Parabeln 2.0. p,,1(%), P?1,.(x) erhält man also die gleiche Güte der Annäherung (in h) Be; 
von x* durch die Parabelnullstellen &,1, &,. wie bei Einschließung von f(x) durch gewisse _ B 
Polynome®) vierten Grades in x. Bei Verwendung von v (x), v(x) hat man anstatt Gleichun- 
gen vierten Grades jedoch nur solche zweiten Grades aufzulösen”). } . 


-. 3. Einschließung der Nullstelle <* von y=[ (&) bei ee der Monotonievoraus- . 
setzung. Es sei L (x) in J(&, <&<2,) eindeutig, f(x) f(®3) < 0, ©* die einzige Nullstelle von N 
: fx) in J,f® (x) (e= 1,2,3) dort stetig. Über die Monotonie von f” (x) sei nichts vorausgesetzt. 2 
Sn Nullstellen der zugehörigen Polynome P1,2(%), P3,1(%) werden also möglicherweise in J z* Rah 
nicht einschließen. Wir werden zeigen, daß sich dann trotzdem mit dem Einschließungsver- ae 
fahren Schranken für x* finden lassen und zwar durch Übergang zu einem passend gewählten 
Gleichungssystem. 
Wir können jedenfalls (dreimal Sri differenzierbare) Funktionen @(x), y(x) angeben, 


Re ara und) 
ist und in 3, <x<x, p’(x) und vw” (x) monoton sind. 


6) a.a.0. }), GI. (24b). 

?) Gleichfalls auf die en. einer quadratischen Gleichung kommt man (siehe R. Beyer, Z. 
. angew. Math. Mech. 16 (1936), S. 345—346), wenn man f (x) durch einen Kegelschnitt annähert, der in einem 
Nachbarpunkt x; von 2* f(x) fünfpunktig berührt. Der Fehler, der damit für z* erzielten N: äherung ist dann 
klein in ö von fünfter Ordnung, wenn man mit ö die Abweichung der Ausgangsnäherung x, von &* bezeichnet. 


PN 912 (2) DRRCHRDEr 
y- a ee 0 0) ER ge ig . 


5: Sp Sp Sur 
s Dabei sei 


Be der Abszisse «; bzw. & -- ie -,) der ER von y—Pa1(®) a var) BE Be! 


Sa. on 0 Tele „are 5) »r „ »,9=P1al2) » Be N = 
i Sı ( ” ” 2] Cr 35 & = — a) ER] \ . LE} Y=931l®) s yayııla) | > 
Sn|» Er > Cr > &r= Sr u 9 2 zu) Er} > Er} y=9,.(2) », y=Y1.[8) 


»* wird (s. Bild 1) von middestens einem der beiden Paare &,, ©; bzw. &,, 211 eingeschlossen 
und zwar, falls f’’(x) in J nicht monoton ist, jedenfalls von %,, 27], wie wir zeigen werden. 

Aus (1) folgt, wenn man darin f(x) durch p(z), y(z) ersetzt und p,,(2) bzw. P,,1(2) 
durch 91,2(@), Yu2(®) bzw. eg Y21(®) 


9,12) —Y12(®) = =9(a)—y(a) ne N Ko—299"6) Ha-ar"n)| 
91,2(2) —Yz1(2) =P(R) — via) — ea) on Er fm, a) 2" (a) + Rz) yP"" (a} | 
P3,1(8) — Y31(8) = PR) Sun ni (X 67 2)? (&;—%) {pP (2) = (3} 


: . (14b). 
P1,2(®) — Y1,2(2) = iR) ul 3, @ — 2) (8. — 2)® {pP (31) u 3} 


In (14a), (14b) sind (bei Beschränkung von x auf das Intervall ©, < x < &,) 31, 3» -.- gewisse . 


Zwischenstellen zwischen &, und &;. Wegen der vorausgesetzten Monotonie von @’ '(2), Y() 
wechseln in unserem Intervalle (x), y"’(x) ihre jeweiligen Vorzeichen nicht: Es besitzen 


‚ also dort 9"’(), y”’ (x) entweder (a) en gleiches Vorzeichen oder (b) durchweg entgegen- 


gesetztes Vorzeichen. (Das bloße Verschwinden einzelner Funktionswerte werde nicht als 
Vorzeichenwechsel angesehen.) Wenn f”’(z) in J nicht monoton ist, kommt (b) nach (13) nicht 
in Frage, so daß der Fall (a) zutrifft. Die Gl. 19 zeigen, daß dann x* durch %, Zr einge- 
schlossen wird. Ist im Falle (a) überdies (9’ (23) — y”’(z,)) durchweg positiv oder durchweg 
negativ für alle 2,, 2, in 2, <2, 23 <%% ..., so folgt aus (14b), daß auch &,, &;r obere und 


untere Schranken für «* geben. Im Falle (b) schließen (nach (14b)) x), 2, die gesuchte Null- 
stelle x* ein. 


‚(Bei p (a) =0 bzw. = (0) BEO ie = P2® Ally = Pa ir bzw. y(a) = ya) = Yya1l®); 


2 ++, Dr DZW. &, 2.2, Er ergeben sich (als sonen EN Schnittpunkte der FinschHe. 
Bungspolynome y = 9,,1(%), y= Ya,1(%), ...) aus den zwischn -A<&E<h gelegenen Null- 


stellen von Polynomen B,,(6) — Fl) =0,...,d,.(8) — Pyı (&)=0 (e Be fer: a 
also als Lösungen von quadratischen Gleichungen der Form 


kot2kulöih) HE, (EI =0 Ne a ee (15). 


(1 1a), pP 


EN 


2 ie 2 Y len 


Br u I Ic: = } = re FA .® ir “ iD er = 
3ER BER Ein rs er u 1 
ES BT TARR ka ie (m — Yin) \ | 

“ £\ | — N 


In (16) ist. om 9 (am): ae a 


Die speziell für y=0 aus ( 13) bis (16) SR Formeln ee nt den träher®) ( 
a den Fall einer monotonen Funktion So ) ns ANSCRnCkEN überein. SEN 


= E Eee Hnschließun vn y=/(@) dureh Parabeln A. der. @. e= en € yr. 


Bm = EI 2a eindeutige Funktion y= f(x) besitze dort genau eine Nullstelle a*, es sei 
Fe) ft (2) = D- Me sei in 9, <x<a; f(r) dreimal stetig nach x differenzierbar und s 


AR END IUTE u EEE 


Es bedeute Yo, = de Meer Yr- Min (@ Y y2) fe —= Max (y4, 9). Wegen 
( =on besitzt y=- le Sr in yı< u Us; N Umkehrfunktion = g(y). Esseien = = Ioıl (Y). 


= Ka n; y bzw. = q,.(y) = B> DR 2.49 die Polynome (höchstens) ‚2ten Grades in nr die auyee 


R i= 
BR durch die Punkte a (21) )), (2 f(xs)) gehen und die Um- 
K. kehrfunktion 2 —=g(y) in (x, Fa) bzw. (&,, f(%,)) in 
1. Ordnung berühren Hi Bild 2). 95,1(%); 91,5(Y) schließen dann 
inyı <y<Yy, die Umkehrfunktion g9(y) ein, wenn. dort 
....%=9(y) dreimal stetig nach y differenzierbar und 


I 


97 . 2 7 
= RN. 7 E MOHGLON ar. > ar LO) 
£ | ay?. Bet 
. ist. Bei unseren ie über f(x) ist die erste 


Bedingung erfüllt. 
Man erhält so &* einschließende Schranken &%,, &ır SO- 2 
fort durch g,,,(0), g1,. (0), zu deren Bestimmung ersichtlich Al 7x, x,) 
keine quadratischen Gleichungen aufzulösen sind. Anderer- Bias 
‚seits ist, wie aus Stetigkeitsgründen für A—0 folgt, nur für 
kleine % gewährleistet, daß x,, %ı ins Innere des Ausgangsintervalles (x}, x,) fallen. 


eh 2 
nel ER ae an 3 a ara or y 


| f# | | < 
RATEN na 
ca Va le re rel 
Mit ©, — a = 0;h, « — 2, =0,h,0< 9,; 9,<2, 0,40, =2 folgt hieraus für kleine A 
up} 1 BES: } ? 
Dyg 0% h? 5 a +: D Ch BR (SR TE 
wobei 0 z0,-0,<H ist: 


8) 2. 2.0. )). 


(a—a-2h. . (19). 


bzw. Se welche g(y) in den Punkten (z1, f(z)) bzw. en 


| re %, eindeutige und viermal stetig nach x differenzierbare Funktion f(x) besitze 


a Bei, 


N das im Ara 2 betrachtete ne fia) = —3: 2 „4, a 
man aus (19) = 14224... << 1,4422490 . er Bi 


Bei-Annäherung von = = 9(Y) dach Taylorsche Polynome NY 
berühren, ergeben sich wie es sein muß, die dreigliedrigen 2 en) ri ade For 


k 5 Kazizh Ei 
LER er... 
a a ü mn u % (20). r% 


Y 
schiiich Formeln für die Abweichungen 2, — 2% 2 — 2* der Näherungen oe auf © 
die Literatur 9) verwiesen. 

Die Paare 2, X; aus (19), (20) schließen #* ein, wenn we Voraussetzungen über Ja. 
insbesondere (17), (18) erfüllt sind. 


s 
RN 


5. Näherungsformeln bei 3 Interpolationsstelen %s Pe: --1 RE, 2); #2. Die in. 


dort genau eine Nullstelle «*. Es sei f(z,) f(z,) < 0. 


A Pr 
a) y= P3,,1,1(8) bzw. y= P,,1,,(2) seien die durch die Forderungen 


b Para) = Prrılı) =S 931,12) = fm) =fm 921,12) =Jı 
ZW. Nr 

Pu12(2ı) =S- P1u12(2m) = fm, Puna(®) =f4 Pın2(2) =f4 
eindeutig festgelegten Polynome 3. Grades in z. “ 

Aus der (bei Markoff!") angegebenen) Darstellung für das Fehlerglied fl) —Pr..a..2,(2) 
bei drei Beer ersieht man, daß die in 2=z, sich durchsetzenden Parabeln 
P2,11(®), P1,1,2(2) in a1 <&< %, die Kurve y=f(z) einschließen, falls dort f'" (=) monoton 

1st. 


Die Gleichungen der FRrächhepiheslolynome lauten (= 2 — m, 2h=2,—1)' 
z EN Er 
Y = P,11(&) = & Ö, Er 3 Y = Pina (£) >= B> Yr Rh 
- 0 v=(0 
mit 
60 — fin Y=In 
hf tm sera Ama la a | 
EN y=b, | j 
a ah Ay amt an | 
Bei hinreichend kleirem h besitzen P,, 11(&) bzw. P,,1,.(£) jeweils genau eine Nullstelle 
&,1,ı bzw. Ein —h<E&<h, welche, falls obige Monotonievoraussetzung hinsichtlich N) 
fa 2) erfüllt ist, ©* einschließen. ‘ 
Weiter gilt. (,=f'(aN) #0, Bu=&ııt Le) j 
of | (e)} & 
BREI BR E Zen 2 as er En * 2 HShLAUT SEHE 7 
821,1 9,12 = 1,17 %n2 271 f 1 „ n nr | 
. 


wobei (wegen |}; 1,2/h| < 1) 


Bere (Zus) | < 04 


b) Bei Annäherung der (unter der Voraussetzung f’(x) +0 in .< ® < x, eindeutigen) & 
Umkehrfunktion x —=g(y) von N = f(%) Es Polynome & = 93,11 (Y) bzw. 2= 9,12 (9%), / 


ist. 


welche N die Punkte.(2,, f(2)), (2m f(&m)); (2 f(x,)) gehen sowie g = g(y) in (2, f(® 
bzw. (2, f(&,)) in 1.0. berühren, erhält a (vgl. hierzu auch Abschn. 4) die a 
r i l Im Fin 
I Sn m er + +-f2 en 
1l er Une BE N (fı — ff: — fm) -\f4+ af Ta er IH, 


I—hrr rat ae 
9) Fr. A. ae 7. angew. Math. Mech. 18 (1938), $. 200, Gl. (35). 
20) 44 .9,.102.2)5 a a 


ee. ur 


Be 


E nr . 
t 5 = SAIEYTM 
" } 1 ’ h N 
AL . ” 2 


wobei 
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KLEINE MITTEILUNGEN Be 


Über die Methode von Graeffe. 


Bekanntlich besteht die Graeffesche Methode zur. 
Lösung algebraischer Gleichungen darin, aus der ge- 


gebenen Gleichung f (x) = 0 mit den Wurzeln ee 


En eine Andere herzuleiten f, (2) =0 mit den Wurzeln 
2, daraus eine neue f, (2) = 0 mit den Wur- 


zeln ©, ..., 2} usw. und diesen Prozeß der Wurzel- 


Br quadrierung so lange fortzusetzen, bis die Gleichung . 


fx (X) = 0 bis auf die gewünschte Genauigkeit in die 


Gleichungen _ 


4,9) =0, 4,&)=0,.:.., H&)=0 


mit den Gradzahlen n,, n,, -: 
+n;=n, zerfällt. Die Wurzeln jeder Gleichung Ms 
haben denselben Modul R;, so daß die Wurzeln sind 


Xst = Rs (cos Ast +isin Ast), t=], ... ne. 


- Daraus ergibt-sich R, sofort als die n;-te Wurzel aus 


dem absoluten Gliede von Ms; = 0. Ebenso ergibt 
sich der Modul r; der zu X5s,: gehörigen Wurzel %s,t 
von f(z) = 0: 

%s,t = Ts (COS Qs,t + Ä Sin az,t) 
als 


nn er k 
s=YR, wo p=2", % 
Alle bisherigen Rechenoperationen verlaufen ganz 
einfach, Unangenehm wird die Rechnung erst für die 


Bestimmung der Argumente as,;. Und hier teilt sich 
die Methode in zwei Modifikationen. 


I. Die erste Modifikation wurde schon 
von Gräeffe selbst angegeben. Er löste nämlich die 
Gleichung M5; (X) =0 uiemals vollständig, sondern 
er bestimmte nur den Modul AR; (und daraus rs). 
Um jedoch die ‚Argumente As, zu finden, löste er die 


Y Gleichung 


en 
d.h. er bestimmte wiederum nur die Moduln r! ihrer 


Wurzeln x’ auf die oben beschriebene Weise. Um die . 


Argumente zu finden, ging er jetzt: folgendermaßen vor. 


Er ordnete, die Moduln rs der Größe nach, dabei 
jeden Modul r; so oft hinschreibend, wie seine Vielfach- 
heit ns beträgt: PR 

01» Q2> ++» On 
Ebenso ordnete er die Moduln r; der Größe nach: 


Ki ’ 
01, 02 +++, On > 


N, wobei nı -H - = ' 


Dann bewies er, 


A, = al &=-lu tb. | 
Dazu muß 6: so Nee sein, daß es höchstens gleich 
der kleinsten Differenz zweier aufeinanderfolgender 
Moduln von f(z) = 0 ist: 

bSmin(r,,,—r) für alle i. 
Denn in diesem Falle schneidet jeder Kreis C'; mit dem 
Mittelpunkt 0 und dem Radiusr; (auf welchem sich 


"ns Wurzeln von f=0 befinden) nur einen einzigen . 


Kreis C7, mit dem Radius &; +Db und dem Mittel- 
punkt —b. 


Ostrowski!) bewies dann, daß es nicht einmal: 


notwendig ist, b genügend klein zu wählen. Im Gegen- 
teil ist b willkürlich. Doch muß man dann die Mo- 
duln oe; von |v; + b| nicht mehr er Größe nach ordnen, 


sondern so, daß 


Eu, Nor: 


(wo c„ das absolute Glied von f (x) ist), eine Beziehung, 
die sich einfach beweisen läßt. 

Aus g: und 0; ergibt sich dann einfach das Argument 
di von &i. 


I. Die zweite Measkatien löst jede. 


Gleichung M5 (X) = 0 vollständig, d.h. man bestimmt 
nicht nur den Modul Rs, sondern auch alle Argumente 
As,ı der Wurzeln X. Dann aber ist bekanntlich das 
Argument @s,t von x nicht, mehr ohne weiters ein- 
deutig, vielmehr hat man es unter 2% Werten zu be- 
stimmen. Wie man die Argumente As,+ auf die ein- 
fachste Weise bestimmt, werden wir unten angeben. 
Zur Bestimmung der as,: aus den As: jedoch hat 
Carvallöo bekanntlich folgende Methode ange- 
gegeben. 

Er bestimmt sukzessive (wenn wir der einfacheren 
Schreibweise halber von jetzt an die Indizes bei den 
X, x, R, r usw. weglassen) . 


LE inc 

VX, YX, YX, ..., VX, 
indem er jeden gefundenen Wert: in die vorhergehende 
Gleichung der Kette einsetzt. Spaltet man somit 


1),Ostrowski: Recherches sur lamöthode de Graeffe et 
les zeros des polynomes et des series de Laurent, Acta Mathe- 
matica, 72, 1940, p.242. In der dortigen Formel (88,1) und der 
vorhergehenden befinden sich zwei Flüchtigkeitsfehler. 


daß für hinreichend kleines b> 0. 
zwei einander entsprechende ‚Moduln gi, 0; zu der- 


selben Wurzel z; von f(z) = 0 gehören, d.h. daß 


ungerader Ordnung: 
ae) +rhe) =, 
so erfüllt x die Gleichung v 


© = — g(a®)/h(aR)., 
Nun ist aber :«? die Wurzel der vorhergehenden 
SE fs (X) = 0, die schon vollständig aufgelöst 
ist. Somit ist 7 


"z=/X=—g(X)/h(R). 
D.h. um die Wurzel der vorgehenden 


die Wurzel der Gleichung £r =0, wel- 
che schon bekannt ist, in die ÄAus- 
drücke g und kh von fi , einsetzen. 
EN Ä Wenn auch g und Ah nur etwa die Gradzahlen n/2 
Ri . haben, so ist doch dassukzessive Einsetzen 
% aller Wurzeln von f£k =0 in alle vorgehenden Glei- 
chungen der Kette sehr mühsam. Daran ändert 
sich auch nichts, wenn man die Rechnungen nur auf 
2—3 Dezimalen oder noch weniger ausführt. 

Damit ist aber die zweite Modifikation der ersten 

a‘ \ in ihrer Zweckmäßigkeit weit unterlegen. Und dies 

©; war-der Grund, daß man den Kunstgriff von Car- 
SB ı '  vallo in seiner bis heute unverändert gebliebenen 
Form vernachlässigte und die erste Modifikation be- 
, vorzugte. 
PN, Nun läßt sich aber der Kunstgriff sehr vereinfachen. 
. x Zunächst wird man überhaupt das Verfahren nicht 
schematisch durchführen, sondern sich erst über- 
zeugen, welehe Moduln r; zu reellen 
Wurzeln gehören. Dies ist schnell festgestellt, 
wobei man sich verschiedener Hilfsmittel bedienen 
kann, und die betreffenden Wurzeln scheiden für das 
Verfahren nach Carvallo aus. 

Bei dem Kunstgriff selber hatte man nun bisher 
nicht beachtet, daß es in Wirklirhkeit überhaupt 
nicht nötigist, irgendeine Division 
auszuführen, ja nicht einmal die 
Funktionswerte g(X) und h(X) zu be- 
rechnen. Man hatte nämlich ganz übersehen, daß 
ja die Wurzel & von f,. _, = 0 nicht völlig unbekannt 


ist, sondern daß man nur zwischen zwei bekannten 


Werten zu wählen hat, nämlich, wenn 


; X=R(cos A-+isin A): 
zwischen 


,,= YR (cos A/2 + isin 4/2) und 2, = —r,. 
Zu einer solch einfachen Entscheidung braucht man 
aber nicht g(X) und A (X) zu berechnen und erst 
recht nicht beide durcheinander zu dividieren. Denn 
das Argument kann nur zwei Werte annehmen, die 
sich um zz voneinander unterscheiden. Man kann bei 


der Bestimmung von © also eine Ungenauigkeit: be- . 


gehen, wenn diese nur absolut genommen kleiner ist 
als . Nun ist 


X) = Here jA+iNeRsinjA= cea 
h(X) = 3,0" Rmcosm A -+i S} c", Rm sin m A =’ 68, 
daher 

C 


= —g(K/hK)=— 7, 


i a8) ? 
6 


d.h. —r hat das Argument &— ß, wo nach obigen 
Ausführungen & nur bis auf 7/2 genau bestimmt zu 
werden braucht und ebenso ß. Dies ist aber sicherlich 
der Fall, wenn man die Quadranten kennt, in welchen 
die Funktionswerte g (X) und h (X) liegen, und dazu 
wieder brauchen wir von obigen: vier Summen nur die 
Vorzeichen zu kennen. Dies ist aber schnell geschehen, 
da die Gradzablen von g und h etwa. n/2 sind und zudem 
in beiden Summen von g(X) jedesmal die Faktoren 
c; Ri die gleichen sind, und ebenso bei h (X). 


ce 


Gleichung f, ,=0 zu finden, hatman. 


ann ist. die Quadr, 


Gm der, wm 


Das Argument: von z ist dann derjenige der beiden 
Werte / bi, 


4/2 oder, at4/2, N 
Angabe in Übereinstimmung ge- 


weicher mit Me 
bracht werden kann, Re: 
Mit dieser Ni von 


der erstengleichwertig. ver 
Es handelt sich also jetzt noch darum, für die 


vollständige Lösung jeder Gleichung 


M(X)=0 von (X) =0ein zweckmäßiges 
Verfahren anzugeben. Dies habe ich in einer 


früheren Arbeit getan?), wovon ich hier nur die R- 


sultate anführe. 
Sei 
M(X)=(0,Xmt+..- +E®=0. 


Dann verfährt man folgendermaßen: 


1. Der Modul R aller Wurzeln von M =0 ergibt 
sich aus ; 


Bm — |Cm/Chl. 


Allgemein ist, wie man leicht beweist, der Quotient 
zweier symmetrischer C absolut, genommen gleich 
einer Potenz von R: > 


RM—2T 2: ICm-r/Er} 3 N 


so daß, wenn m ungrade ist, R sich einfach als Quo- 
tient zweier mittlerer Glieder ergibt: 
R=](0,,,/Cl, wo m=2s+1. 

2. Nachdem man somit, R bestimmt hat, setze man 
X=RY, wo Y=e” unbekannt ist. Setzt man 
diesen Ausdruck von X in M(X)=0 ein, so erhält 
man, worauf man bisher noch nie geachtet hatte, 
eine reziproke Gleichung M"(Y=-0inY, 
weil mit e'# auch eF _ y-ı eine Wurzel ist. 

3. Man dividiert M’(Y) so oft durch Y—1 wie 
dies möglich ist, etwa m’mal. Dann sind m’ Wurzeln 
von M(X)=0 gleich R. } 

4. Der Quotient Q(Y) = M’/(Y—1)”" ist ebenfalls 
ein reziprokes Polynom. 

Setzt man daher Y + Y’l=Z=2cosW in Q(Y) 
ein, so bekommt man eine Gleichung 9’ (Z) = 0 vom 
Grade m’’ = (m—m')/2, deren Wurzeln alle reell und 
kleiner als 2 sind (absolut genommen). Diese Gleichung 
läßt sich also einfach nach ‚Grae£ffe lösen und möge 
die Wurzeln Z,, Zu +++, Zmr liefern. 

5. cos Wi = Zi/2rliefert dann m’’ Werte für W, 
also 2m’’ Werte von Y. 

6. Diese 2m’’ Werte liefern zusammen mit derım’- 
fachen Wurzel X = R das vollständige System X=RY 
der Wurzeln von M (X) =0. 

Diese Methode der Lösung von M (X) = 0 und im 
Anschluß daran von f(«)=0 versagt in kein- 
nem noch so komplizierten Falle. 

‚Man hat zahlreiche Kunstgriffe vorgeschlagen, um 
die Konvergenz des Graeffeschen Ver- 


N ®) B. Bodewig: On Graeffe’s Method for Solving Algebraic 
EAN Quaterly Applied Matematics, USA,, vol.IV, 1946, 
P.177—190. 


| 
| 
Carvallo ist die zweite Modifikation dr; 


leichungen M,=0, ..., M;=0 


ezug auf die Wurzeln x Ban tökineh. 

r wäre es sehr unzweckmäßig, die geleistete _ 
eiseite zu schieben und erst f(@a)=0 zu 
ransformieren und auf die transformierte Gleichung die 


Graeffe-Transformationen anzuwenden. Denn man 


um einige transformierten Gleichungen zu verlängern, 
spalten zu sehen. 


Wie ich nämlich in genannter Arbeit: bewiesen habe, 


©. konvergiert die Folge der Werte 
8 = r Fa Ei VRs. ME 


EN tür Bl, 2... 
gegen rs. Dabei verstehe ich unter quadratischer 


“wert X, daß 
ba % ' a $ Al 
rue Zusammenfassung. 

1. Dieschematische Ausführung der Graeffe-Trans- 
formationen bis zur Aufspaltung der transformierten 
Gleichung ist allen Kunstgriffen überlegen, da diese 
Methode quadratisch-konvergente Folgen liefert. 
IR 2. Der oben gezeigte Weg der Lösung jeder Glei- 

chung M;(X)=0 gefolgt von dem vereinfachten 
 Kunstgriff von Carvallo ist die einzig sichere, 
' nie versagende und dabei vorteilhafteste Methode, 


Den Haag. - Bodewig. 
- Beitrag zur inversen Interpolation. 
- Wir gehen aus von den Identitäten!) 


== 204 (Yo) [ro 21] 
— Yy—Y) (YıY) lRorız2l (1) 
+W— 9) yı 9) We Y)larorı22]) 
r un 
E 2=2+ (y— Yo) [&o 2ı] 
>> 9) (yı—Yy)le-ı Fo 83] (2). 
4 — (y— 4-1) (y — 90) Yı—W)ler-ırori] 
A Hier ist y_4 < Yo. <Y <Yı < Y, angenommen. Wegen 
1 (dx 
} raal=z Ay)’ Yo <Nı <S Ye 
| _ 1 (die | 
B#. | [r2_1 or] = 31 dys/m’ Y_-ı < Na <Yı 
ö erkennt man, daß die aus den Gleichungen (1) und (2) 
folgenden Approximationspolynome 2. Grades den 
j gesuchten Wert einschließen (vorausgesetzt, daß der 
-dritte Differentialquotient im Intervall y_,, y, das 
Vorzeichen nicht wechselt). Ausrechnung und Um- 
formung ergibt dann für die inverse quadratische 
Interpolation die Einschließungsformeln ?) 
1) Hier bedeuten 
j ; '% 7% J 
[2 %] >= Bo usw. 
. die ‚„‚Dividierten Differenzen‘. In meiner Schreibweise ist also 
Me [vor] = ® = usw. 
im Gegensatz zu der sonst üblichen. 
’ 2) Diese Formeln wurden vom Verfasser bereits auf anderem 


wege abgeleitet (vgl. diese Zeitschrift Bd.28 (1948), 3.122). 


Bet 
? 


mäßig. Denn bevor man | 

ber die Wurzeln hat, hat man — 
bgesehen— mehrere Graeffe- | 
ühren müssen und hat erst |v—=x, 
Aufspaltung der Gleichungen | 


. brauchte nur die abgebrochene Gleichungskette noch _ 
' um die neue Gleichung ‚fr (X) =0 sich völlig auf- 


und festes s quadratisch. 


© Konvergenz einer Folge x,, &;, ... mit, dem. Grenz- 


. ergeben Formel (5) und (6) 


ya nal nn 


Hier bedeutet & Le 
Für die Anwendung i 


= 2 AyAyı  . 


EL A| Bea 
Br ER Er -2 Ayıdyo PP, 


und erhalten so die einfacheren Einschließungsforme 


‚Als Beispiel soll 
u EEE } 

i gran v2 (4 Y2= 0,70710 oia12) N 
berechnet werden. Mit den Funktionswerten 


sin 0,784 = 0,70611 74395 sin 0,786 — 0,70753 22158 
sin 0,785 = 0,70682 51811 sin 0,787 — 0,70823 85490 


.0,78539 81633 < = < 0,78539 81636 . 


Der Mittelwert ist bis zur 10. Dezimalstelle richtig.” 
Die Formel für den Mittelwert von (5) und (6) lautet 
mm nn 


Yyı—y Ay au m. 
4 Ay, Ayı 


' Durch Anwendung der Näherungsgleichung 
Y— Yo) Ayı Sue 


= Et ® 
[4 


Y—Yo h_ 
Ayo \ 


W—y)Ayı + mı—y)Ay Ayo BR 
We) nı-y) An—Ayı) TR 

Ay, AyoAyı r ER 

können wir Formel (5) schreiben i x je 


(yo) Ayı ‘a ER 
(y—Yy0) Ayı + Yı —Y) Ayo i Ir j 
y— Yo Yı 4 Ayo — AYı 
Ayo 2° ‚Ayodyı 


d.h. die inverse lineare Interpolationsformel 


et r@ 


+0 


R (4 — Yo) Ayı Eu 
(y— y) Ayı + (yı—y) Ayo 


(8) 


s=%+ 


ergibt (bis auf Glieder 3. Ordnung) einen entgegen- 
gesetzt gleichen Fehler wie die Regula falsi 
UN 0 
=1,+ 73, 
und eignet sich daher vorzüglich zur Fehlerabgrenzung 
der inversen linearen Interpolation. 


Vertauschen wir in Gl. (1) x, und y, mit , und Yo, 


so folgt 


Fa — yo) [leo 81] 
: nina | 9) 
— (y— yo)? (yı —Y) [ex 20 &ı] en 
und wegen & z 


1 dx. 
a ee I 1-(5),) 
ergibt sich das Approximationspolynom 2. Grades 


"3 ae 1 
rar yzylsmie, Hin) = (10). 


« .Y177%90 Yı Yo 


Vertauschen wir dagegen in Gl. (1) x, und y, mit z, 
und y,, so folgt 
2 = 0 + (y— Yo) [@otıl 
un Nlenzlı (1) 
+(y— Y0) Yı— Y) Leo 21 81) 


oder, wegen « 


Nee Fig 
Haag, (5), ee 


ERRBERE Kst ae BE RAR 7278 
3 | Pe a a rr] 


Aus (9) und (11) erkennt man, daß die Approxima- 
tionspolynome (10) und (12) den gesuchten Wert 
einschließen. Für ihren Mittelwert erhalten wir die 
nicht nur genauere, sondern auch einfachere Formel 


N a I RAR zu (3). 


Ihre, Genauigkeit wird bei & = 10”? im allgemeinen 
bis zur 10. Dezimalen reichen, 
Als Beispiel soll wieder 


T =aresin = y2 (4 Y2 = 0,70710 o7s12) 

berechnet werden. Mit den Tafelwerten 

sin 0,785 = 0,70682 51811 cos 0,785 = 0,70738 82692 

sin 0,786 = 0,70753 22158 cos 0,786 = 0,70668 10904 

ergibt: unsere Formel bis zur letzten Stelle richtig 
m 


7 = 18539 81634. 


BUCHBESPRECHUNGEN 


Dr. Wolfgang Haack, Differentialgeo- 
metrie. Teill. (Bücher der Mathematik und Na- 
turwissenschaften, herausgegeben von Dr. Henry 
Poltz.) 1368. m. 19 Abb. Wolfenbüttel und Han- 
nover 1948. Wolfenbütteler Verlagsanstalt G.m.b. H. 
(Notdruck) Preis brosch. 9,— DM. 

Nach einer sehr knappen, vom Invariantenbegriff 
ausgehenden Einführung in die Vektoralgebra behan- 
delt der Verfasser die Raumkurven, die Elemente der 
Flächentheorie (Fundamentalformen, Sätze von Meus- 
nier und Euler, Gauß’ theorema egregium), die Abbil- 
dung zweier Flächen auf einander (längen--, winkel- und 
flächentreue Abbildung), die Ableitungsformeln von 
Gauß und Weingarten und die Verträglichkeitsbedin- 
gungen von Mainardi und Codazzi, die Geomterie auf 
der Fläche (geodätische Linien, Parallelismus von Levi- 
Civita, Gauß-Bonnetsche Integralformel, Poincares 
Modell der hyperbolischen Ebene) und schließlich die 
Minimalflächen. Die Fülle des auf 133 Seiten zusam- 
mengedrängten Stoffes bedingt eine erhebliche Kürze 
der Darstellung, die dem Anfänger vielleicht zu schaf- - 


IT se t 
Pr hie cr “oe h = 
>| "ee, ya 

KO-l—n) ytant+zor dw n- 
abgeleitet. Mit 5 7 RE 
Mr wu 


’ + 
AR 3 
0,7% ; Ar 


können wir dieser Gleichung die Form geben. | 
EM 1 } 5 ee 
y=9Yyat = Pi Inn ee 2 


4 an | 
und hieraus folgt nach Vertauschung von zmitysofort 
die re für die inverse quadratische Inter- 

lation. er 
TB Be it sei auch darauf aufmerksam 
ht, daß die vom Verfasser früher angegebene 
Formel für die inverse lineare Interpolation dee 
Exponentialfunktion ®) 
= +yet%—]) R 
das erste Glied ist einer allgemeinen Entwicklung 
nach Potenzen des Klammerausdrucks. Das ergibt 
sich sofort aus der Potenzreihe 


oo 8 
hi+)=-I- tt, 
g=1 


8 
wenn man x durch 


2 
Yo 


ersetzt. Für y= e”, yy = e”* entsteht dann die Ent- 
wicklung 


& =a+Wwen—)— 5 yernr—1l)?+... (14). 
Berücksichtigt man nur die beiden ersten Glieder, so 


erhält man für die inverse quadratische Interpolation 
der Exponentialfunktion die Formel 


z=0+ 4 werr—1){2— yer—1} (15). 
Der Fehler ist kleiner als 
wenn. 
F.K. Rubbert. 


®) Formel (2,3) in der o0.a. Arbeit des Verfassers, die dort aber 1 
einen Druckfehler enthält, indem es nicht 


Göttingen. 


e-?, sondern &-% 
heißen muß. Entsprechend ist (2,4) zu berichtigen. 


fen machen wird, obwohl sich das didaktische Geschick 
des Verfassers auch hier wieder bestens bewährt. Er | 
knüpft seine Betrachtungen, wieesindeneinführenden ' 
Vorlesungen über Differentialgeometrie üblich ist, an 

die anschaulichen geometrischen Sachverhalte an, ver- | 
steht es aber trotzdem, den Leser schon an den invari- 
antentheoretischen Kern der Probleme heranzuführen. 
Erfreulich sind die zahlreichen historischen Anmerkun- 

gen und die Ausführlichkeit, mit der der Verfasser im 
Hinblick auf die Anwendungen in der Kartographie die 
Abbildungen der Flächen aufeinander behandelt. Miß- 

fallen hat dem Rezensenten nur, daß der Verfasser dem 
Einbruch des Zeitungsjargons in die wissenschaftliche 
Literatur Vorschub leistet, indem er von „‚Bertrand-- 
Kurven“ und „Darboux-Vektor‘‘ spricht und eine 
Funktion ‚‚nach Taylor‘ entwickelt. Bei seinen wei- 

teren Arbeiten möge der Verfasser beherzigen, was in 
Wustmanns Sprachdummheiten über dieses Thema ge- 

sagt wird! Dann sei noch bemerkt, daß man bei der 
Taylorschen Entwicklung eines Vektors nicht von 
einem Zwischenwert din der Restgliedformel spre- 


et ei £ ir ‚ Mer 7 = 3 5 
rechnung in Mathematik, Physik und Technik zu 
‚behandeln. EN BEN... 

Das Bändchen gibt eine gute Einführung indas 
behandelte Gebiet. Besonders gut wird es neben Vor 
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ur selbe ist wie die Verfolgungskurve — ‚„‚Hundekurve‘“ lesungen zu gebrauchen sein. Für das Selbststudium ° 
Kr nennt sie der Verfasser! Das Buch zeigt bei aller zeit- scheint mir die Darstellung an manchen Stellen Ö: 
® bedingten Schlichtheit der Aufmachung doch ein vor- etwas zu konzentriert zu sein. - N 
' bildlich klares Satzbild und vorzügliche Figuren. Alles Dresden. j “Willer so 028 


Freiberg. az Grüss. i Karl Daeves und August Becke], Großzahl- 
EN s Rh orschung und Häufigkeits-Analyse 
Sn  nurIng. habil. H. Athen, Vektorrechnung Ein ragen, 65 8. m. 17 Abb. Wen 


(Bücher der Mathematik und Naturwissenschaften, _ Bergstr. und Berlin 1948. Verlag Chemie GmbH. 


herausgegeben von Dr. Henry Poltz2). 908. Preis brosch. 3.20 DM. a 
i 3, r Et 
m. > Abb. Wolfenbüttel-Hannover 1948. Wolfen- Das Bändchen gibt eine Anleitung zur Durch- 
ag nal GmbH. (Notdruck.) Preis führung von Großzahluntersuchungen insbesondere 
de: nr DM, PATE RR unter Verwendung von Häufigkeitspapier und Wahr- it 
as Bändchen will eine Einführung in die Vek- scheinlichkeitsnetzen, wobei meist logarithmische _ L 


torrechnung geben. In der Darstellung schließt eg Merkmalsskalen benutzt werden. Dabei handelt es. 


sich an die bekannten größeren Werke an, benutzt 
. dabei aber die ältere Schreibweise. (Für skalares 
Produkt runde, für Vektorprodukt eckige Klammern.) 
Nach Behandlung der Vektoralgebra, der Vektor- 
‚analysis und ihrer Anwendung in der Potentialtheorie 
folgt eine kurze anschäuliche Einführung in die 
Tensorrechnung. Nur wenige Beispiele werden ge- 
bracht; denn Verfasser beabsichtigt in einem Ergän- 


sich darum, Kollektive beliebiger Verteilung in eine 
kleine Zahl von Teilkollektiven mit Gaußscher Ver- 
teilung aufzuspalten. Diese werden durch Mittelwert 
und durch die Spanne der Klassengrenzen charak- _ 
terisiert, außerhalb deren beiderseits je 5% der 
Kollektivgegenstände liegen, also nicht wie üblich 
durch Mittelwert und Streuung. de 


Dresden. Willers. 


 zungsband die wichtigsten Anwendungen der Vektor- 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen: 
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"Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleib, vorbehalten). er 


sik der Kolloide. XVI+2378. mit 1 Titl- ° 
bild, 143 Abb. und 34 Tab. Dresden und Leipzig 1948. E 


: 
; 
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‚Dr.-Ing. Fr. BoSnjakovi6 (o. Prof. an der Universität 
Zagreb), Technische Thermo-Dynamik 


(Wärmelehre und Wärmewirtschaft in Einzeldarstel- 

lung. Herausgegeben von Dr.-Ing. W. Pauer, Bd.11, 
I. Teil) 3. unveränderte Auflage, XVI + 327 8. mit 
249 Abb. und 3 Tafeln. Dresden und Leipzig 1948. 
Verlag Theodor Steinkopff. "Preis: geb. 18,—; brosch. 
17,— DM. 


Dr.-Ing. Hans Umstätter (Kaiser Wilhelm-Institut 
für Physikalische Chemie und Elektrochemie), Struk- 
turmechanik. Ein Beitragızur Phy- 


Verlag Theodor Steinkopff. Preis: gebunden 19,—;. 
broch. 17,— DM. . 


Dr.-Ing. Johannes Johannson, Das Cross- 
Verfahren, Die Berechnung biege- 
fester Tragwerke nach der Methode 
des Momentenausgleichs. 123 S. mit 18 
Zahlenbeispielen und 137 Abbildungen. Berlin— 
Göttingen— Heidelberg 1948. Springer-Verlag. Preis: 
brosch. 14,40 DM. 


- NACHRICHTEN 


Richard Grammel sechzig Jahre alt. 


Esist eingroßer Kreis von Menschen, die des Jubilars 
freudig bewegt gedenken: Die Fachgenossen in aller 
Welterkennen den hervorragenden Gelehrten und For- 
scher schon seit langem als einen der ersten in ihrem 
Kreise an, eine große Anzahl von Schülern bewundert 
und verehrtihren trefflichen Lehrer, die Kollegen und 
akademischen Mitbürger danken dem tatkräftigen Mann, 
der ihre Hochschule aus Trümmern neu begründete, 
seine Freunde lieben den warmherzigen Menschen und 
schätzen den geistreichen Erzähler und Plauderer. 

Richard Grammel hat eine bedeutende Wirksam- 
keit entfaltet; sie reicht weit über den Kreis seiner 
eigentlichen Schüler hinaus. Es gibt wohl kaum einen 
Ingenieur, und gewiß keinen wissenschaftlich inter- 
essierten, dem der Name Grammel nicht ein fester Be- 
griff wäre. Dieser Nameistinsbesondere verbunden mit 
drei großen Publikationswerken: Mit der Zeitschrift 
„„Ingenieur-Archiv‘, die Grammel gründete, und die 
er seit nahezu zwanzig Jahren leitet, mit dem ‚‚Bie- 
zeno-Grammel“, der zur ‚‚Bibel‘‘ des Faches ge- 
worden ist, und— etwas weniger in die Breite wirkend 
— mit der Redaktion der drei Mechanikbände des 
„blauen‘‘ Handbuches der” Physik. 

Im, weiteren wissenschaftlichen Werk des Jubilars 
beeindruckt den Betrachter die Konsequenz, mit der 
Grammel die Probleme eines Gebietes aufgreift, zu 


Ende denkt und abschließt, wie er\dann zu einem Nach- 
bargebiet übergeht und dort dieselbe gründliche und 
vollständige Arbeit leistet. Die etwa sechzig Schriften, 
die das Verzeichnis nennt, lassen sich nahezu vollständig 
in einige wenige Gruppen einordnen: Nach einigen Ver- 
öffentlichungen aus verschiedenen Themenkreisen 
wandte Grammelsich dem Kreise | zu. Es entstanden 
eine Reihe von wesentlichen Untersuchungen, die ihre 
Krönung fanden in der Monographie ‚‚Der Kreisel“ 
vom Jahre 1920. Die Untersuchung über ‚Die kriti- 
schen Zustände rasch umlaufender Wellen“ bildete den 
Übergang vomKreiselzu den Stabilitätsproblemen, 
die die Jahre bis etwa 1929 ausfüllen. Dann leitete die 
Arbeit“,Über den Einfluß der Wellentorsion auf die 
kritische Drehzahl“ über zu der großen Gruppe von 
Arbeiten aus den dreißiger Jahren über die Torsions- 
schwingungen der Wellen von Maschinenanlagen, 
durch die die praktisch so brauchbaren Rechenmethoden 
geschaffen und in eine handliche Form gebracht wurden, 
und in denen überdies eine Reihe grundsätzlich neuer 
Erkenntnisse zutage gefördert wurden, wie sie sich ins- 
besondere um die Stichworte ‚‚Torsion erster Art‘ und 
„Torsion zweiter Art“ ranken. Es schließt sich dann 
die Gruppe der Arbeiten über die Scheiben an. 
Zusammengefaßt werden die Problemkreise, denen 
diese Untersuchungen angehören, in dem, schon ge- 
nannten, gemeinsam mit C.B. Biezeno ‚verfaßten 


Hr 
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uns, so dürfen wir, ohne enttäuscht zu werden, hoffen, 


Werk „‚Technische D 
shr als tausend Sei 


Tatsache neu aufgefaßt undyoriginell dargestellt ist, 
Diese Besonderheit Grammelscher Arbeitsweise, eine 


Frage neu zu stellen oder eine alte auf überrasche 


neue Art zu beantworten, spiegelt sich ferner deutlich _ 
wieder in mehreren Arbeiten, von denen. als kennzeich- 


nend vor allem‘ aus neuerer Zeit besonders 


seien: „‚Kinetisch unbestimmte ER Br- 


klär der hohen Sprengfestigkeit umlau: 
Scheiben“, „Eine Verallgemeine der Kreis- und 
Hyperbelfunktionen‘“ und (die noch im Druck befind- 
liche) ‚‚Scherprobleme’’, j 

Es ist noch lange nicht Zeit, Rückschau zu halten und 
eine Würdigung im ganzen zu geben. Die auf ihrer 
vollen, Höhe stehende Schaffenskraft des Jubilars wird 


noch viele weitere Ergebnisse seines Forschens bringen 
und uns noch wesentliche Erkenntnisse und Einsich 
verschaffen. 


Seinen Lebensweg begann Richard Grammel im 
Schwabenland. Er wurde am 3. März 1889 in Kloster- 
reichenbach im Oberamt (Kreis) Freudenstadt ge- 
boren. Seine Studien in Stuttgart, München und 
Tübingen schloß er im Jahre 1913 mit der Promotion 
in Tübingen ab. Dann habilitierte er sich 1915 in, 
Danzig, wo er einige Jahre Assistent beiHans Lorenz 
war. Aus einer kurzen Tätigkeit an der Universität 
Halle berief ihn die Technische Hochschule seines 
Heimatlandes im April 1920 zum Ordinarius für Tech- 
nische‘ Mechanik und Technische Wärmelehre. In 
Stuttgart lebt und lehrt Grammel seitdem. . 

Die Erfolge des Forschers, Lehrers und tatkräftigen 
Mannes haben eine Reihe von Anerkennungen und 
Ehrungen gebracht. Den erst Vierzigjährigen erkor sich 
seine Hochschule zum Rektor im Jahreihrer Zentenar- 
feier. Die Eidgenössische Technische 
Zürich reihte ihn zugleich unter ihre Ehrendoktoren 
ein. Nach dem Zusammenbruch des Jahres 1945 berief 
ihn seine Hochschule erneut, und zwar drei Jahre 
lang, zum Rektor. Hier führte er den Neuaufbau in so 


hervorragender Weise durch, daß die Technische Hoch- 


schule Stuttgart einen beachtlichen Vorspr vor 
ihren Schwesteranstalten erlangte. Man gab dem Dank 
für diese Leistung Ausdruck durch die Ernennung zum, 
Ehrenbürger der Hochschule. In der Gesellschaft für 


ZUSCHRIFT AN DEN HERAUSGEBER 


Zu G, Plato: Über den Zusammenhang von 
Schwingungsdauer und Rückstellkraft. Z. 
angew. Math. Mech. 28 (1948), S. 94. 


Ein interessantes Problem ist in dieser Kl. M. an- 
geschnitten. „Gefühlsmäßig‘‘ bin ich mit dem all- 
gemeinen Ergebnis nicht einverstanden. Nimmt 
man bei einem schwingenden Einmassensysteme an, 
daß im ' ersten Intervalle des Maximalausschlages 
0 <a<a, die Schwingungsdauer 7, (a) veränderlich 
ist, 60 kann für das zweite Intervall a, <a <a; sehr 
wohl eine konstante Schwingungszeit 7, (a) bestehen. 
Für das erste Intervall ist das Potential dann keine 
Funktion von a von der Gestalt ca? und für das zweite 
Intervall auch nicht,, Für dieses Intervall ist das 
Potential durch den Verlauf im ersten Intervall und 
durch 7, (a) bestimmt. 

Das anschaulich gewonnene Ergebnis werde ich in 
einer getrennten Kleinen Mitteilung beweisen, 


Hochschule in, 
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Göttingen: Zu ordentlichen Mitgliedern « 
demie der Wissenschaftenin i wurden 
Dr. W. Magnus, 'o. Prof. der an 
Universität Göttingen, und Dr. F. Rellich, ebenfalls 
0.Prof. d. Mathematik an.der Universität i 

Dr. Stellmacher hat sich im Herbst 1948 an 
Se Universität Göttingen für das Fach Mathematik 

abilitiert. r . Ber 


Ich überlasse es Herrn Plato, den Fehler in seiner 
mathematischen Beweiskette zu finden; zum Schlusse 
kann doch nur einer von uns recht haben. 


Schlewecke über Bockenem, G. Weber. 


Erwiderung. 

Die mitgeteilte Überlegung beschränkt sich auf den 
Fall eines durch eine Potenzreihe darstellbaren Po. 
tentials, was der Kürze halber zwar nicht ausdrücklich 
erwähnt aber doch aus dem Text ersichtlich war, da 
man ja andernfalls nicht von analytischer Fortsetz 
sprechen kann. Die Übertragung der Überlegung a 
nicht-analytische Potentiale scheitert daran, daß zu 
einer Folge Yn(x) eine Folge Tn (a) nr =1, 2, .. BR 
gehören kann, deren une in einem gegebenen 
a-Intervall gegen Null strebt. 


Darmstadt. G. Plato; 


